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We analyse the sensitivity ol all experimentally observable asymmetries and energy distributions 
for the neutron /?~-decay with a polarised neutron and unpolarised decay proton and electron 
and the lifetime of the neutron to contributions of order 10^'* of interactions beyond the Standard 
model (SM). Since the asymmetries and energy distributions are expressed in terms of the correlation 
coefRcients of the neutron /3~-decay, in order to obtain a theoretical background for the analysis of 
contributions beyond the SM we revise the calculation of the correlation coefficients within the SM. 
We take into account a complete set of contributions, induced to next-to-leading order in the large 
proton mass expansion by the "weak magnetism" and the proton recoil, and the radiative corrections 
of order (q/vt), calculated to leading order in the large proton mass expansion. We confirm the 

■ results, obtained in literature. The contributions of interactions beyond the SM we analyse in the 
linear approximation with respect to the Herczeg phenomenological coupling constants, introduced 
at the hadronic level. Such an approximation is good enough for the analysis of contributions of 

■ order 10~* of interactions beyond the SM. We show that in such an approximation the correlation 
coefficients depend only on the axial coupling constant, which absorbs the contributions of the 
Herczeg left-left and left -right lepton-nucleon current -current interactions (vector and axial-vector 

i-^ , interactions beyond the SM), and the Herczeg scalar and tensor coupling constants. In the lifetime of 

the neutron in addition to the axial coupling constant the contributions of the Herczeg left-left and 
left-right lepton-nucleon current-current interactions (vector and axial-vector interactions beyond 
, the SM) are absorbed by the Cabibbo-Kobayashi-Maskawa (CKM) matrix element. 

PACS numbers: 12.15.Ff, 13.15.+g, 23.40.Bw, 26.65.+t 
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^ ! I. INTRODUCTION 

m 
o 

• . In this paper we propose a consistent analysis of the sensitivity of all observable asymmetries and energy distributions 
^ ' of the neutron /3~-decay n p + e~ + with a polarised neutron and unpolarised decay proton and electron and the 
, lifetime of the neutron to contributions of order 10^^ of interactions beyond the Standard model (SM), described 
at the phenomenological level. Such an order of corrections beyond the SM has been pointed out by Ramsey-Musolf 
and Su 0] within the Minimal Supersymmetric extension of the SM (MSSM). As has been shown in [3|~0i at the 
phenomenological level the neutron /3~ -decay may be described by eight complex coupling constants determining the 
strength of interactions beyond the SM. As a result possible deviations from the SM, causing the contributions of 
^ ' order 10~^ to the correlation coefficients of the neutron /3~-decay and the lifetime of the neutron, may be determined 
in terms of vector, axial-vector, scalar and tensor Icpton-baryon weak coupling constants @. 

It is well-known [ll,[2| that in the non-relativistic approximation to leading order in the large proton mass expansion, 
which is equivalent to the leading order of the heavy-baryon approximation, and in the rest frame of the neutron the 
SM with weak V — A interactions @ describes the /3~ -decay of the neutron in terms of two coupling constants Gv 
and Ga [11,13 (see also [l^). The coupling constant Gy is defined by the product Gy = GpVud of the Fermi coupling 
constant Gp = Qw/^^w 0' where gw and M\y are the electroweak coupling constant and the 14^-boson mass, 
and Vud is the matrix element of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix Q. The coupling 
constant Ga is equal to Ga = XGy, where A = Ga/Gv is the axial coupling constant, induced by renormalisation 
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of the axial hadronic current by strong low-energy interactions For the weak interactions, invariant under time 
reversal, the coupling constant A is real. 

The observables of the neutron /3~ -decay with unpolarised particles are the lifetime of the neutron t„ and the 
correlation coefficient ao, describing correlations between 3-momenta of the electron and antineutrino to leading 
order in the large proton mass expansion in the rest frame of the neutron. Experimentally the correlation coefficient 
ao can be determined, for example, by measuring either the electron-proton energy distribution a{Ee,Tp), where 
and Tp are the total and kinetic energies of the electron and proton, respectively, or the proton-energy spectrum 
a{Tp). As a function of the electron energy Ee the correlation coefficient of correlations between the electron and 
antineutrino 3-momenta we denote as a{Ee). 

The neutron /3~ -decay with a polarised neutron and unpolarised decay proton and electron is characterised also 
by two additional observable correlation coefficients Aq and Bq, describing to leading order in the large proton 
mass expansion correlations between the neutron spin and 3-momenta of the electron and antineutrino, respectively. 
As functions of the electron energy Eg the correlation coefficients of correlations between the neutron spin and 3- 
momenta of the electron and antineutrino we denote as A{Ee) and B{Ee), respectively. The lifetime of the neutron 
and the correlation coefficients under consideration are of order t„ ~ 880 s, qq ^ a{Eg) ^ Aq ^ A{Ee) ~ —0.1 and 
Bq ^ B{Ee) ^ 1 The couplingconstants Gy and A define the main contributions to the lifetime of the neutron 
and the correlation coefficients (see also [l^)- 

However for the description of the neutron /3~ -decay at the modern level of experimental accuracies when the 
experimental data on the axial coupling constant A = —1.2750(9) [l| and the lifetime of the neutron Tn^^^ = 878.5(8) s 
|13 | have been obtained with accuracies 0.07% and 0.09%, respectively, two coupling constants Gy and A are not 
enough for the correct description of the properties of the neutron /3~-decay. Indeed, it is well-known that the radiative 
corrections to the lifetime of the neutron 13l|-j30j. calculated to leading order in the large proton mass expansion 
within the SM with the V — A weak interactions and Quantum Electrodynamics (QED), give the contribution of about 
3.9% [13 (see also [l|). This allows to describe correctly the lifetime of the neutron [l^l. The radiative corrections 
are very important also for the correct determination of the Fermi coupling constant Gp- It may be extracted from 
the experimental data on the weak coupling constant G^ of the /z~-decay jjT — > + + i^e with the account 
for the radiative corrections [l3|-[lBl- The contributions of the radiative corrections to the ratios of the correlation 
coefficients a{Ee)/ao and A{Ef,)/AQ are by order of magnitude smaller compared with the contribution to the lifetime 
of the neutron, whereas the correlation coefficient B{Ee) has no radiative corrections to order a/7r [3llf3^. 

Another type of corrections, which may be calculated within the SM and should be taken into account for the 
description of the neutron /3~-decay on the same footing as the radiative corrections, are the contributions of the 
"weak magnetism" and the proton recoil [S^IIJ] (see also [l^), calculated to next-to-leading order in the large proton 
mass expansion. 

The radiative corrections and the corrections from the "weak magnetism" and the proton recoil, calculated to 
leading order and to next-to-leading order in the large proton mass expansion, respectively, define a complete set 
of corrections to the observables of the neutron /3~-decay, which should be taken into account within the SM as a 
background for the analysis of contributions of order 10"'' of interactions beyond the SM. 

Since the observable asymmetries and energy distributions are defined in terms of the correlation coefficients [l| , for 
the aim of the paper to analyse the sensitivity of the observables of the neutron /3~ -decay to contributions of order 
10~^ of interactions beyond the SM we revise the calculation of the correlation coefficients of the neutron /3~-decay 
with a polarised neutron and unpolarised decay proton and electron and the lifetime of the neutron within the SM 
with V ~ A weak interactions. In addition to the non-relativistic terms ao, ^o and i?o, calculated to leading order 
in the large proton mass expansion, we take into account 1) the contributions of the "weak magnetism" and the 
proton recoil to next -to-leading order in the large M expansion, which provide a complete set of corrections of order 
1/M, where M is an average mass of the neutron and proton M = (m„ + TOp)/2, and 2) the radiative corrections of 
order o/tt, calculated to leading order in the large M expansion, where a = 1/137.036 is the fine-structure constant 
0. The small parameter of the large M expansion is Eq/M ~ 10~^, where E^ ^ 1 MeV is the end-point energy of 
the electron-energy spectrum. The parameter Eq/M ^ 10~^ is commensurable with the parameter of the radiative 
corrections a/ir ^ 10~^. As we show below due to the strong dependence on the axial coupling constant A the 
numerical values of the 1/M corrections vary from 10^^ to 10^^. In turn the contributions of the radiative corrections 
do not depend on the axial coupling constant A. Following Sirlin [l8| we show that the contributions of the radiative 
corrections, depending on the axial coupling constant A, may be absorbed by renormalisation of the Fermi Gp and 
axial A coupling constants. The contributions of radiative corrections to the ratios a{Ee)/aQ and A{Ei,)/Aq are of 
order 10~^. The correlation coefficients, calculated within the SM to order 1/M and a/n, determine the theoretical 
background for the analysis of contributions of order 10"** beyond the SM. 

A phenomenological analysis of contributions of interactions beyond the SM model shows that these interactions 
induce only 1) the energy independent contributions and 2) the contributions proportional to me/£'e, where me is the 
electron mass. We show below that the contributions of the "weak magnetism" and the proton recoil into the terms 
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proportional to nie/Ee of the correlation coefficients a{Ee), A{Ee) and B{Ee) are of order IQ-^ - 10"''. Thus, the 
obtained theoretical expressions for the correlation coefficients should be taken into account for a correct experimental 
determination of contributions of order 10~* beyond the SM. 

The radiative corrections to the correlation coefficients a(i?e)/ao and A{Ee)/AQ determine the most important and 
complicated part of the corrections of order 10"'^. In the neutron /3^-decay the radiative corrections are defined by 1) 
the contributions to the continuum-state /3~-decay mode from one-virtual photon exchanges, VF~boson and Z-boson 
exchanges and QCD corrections 124 . l27l- (30l| and 2] the contribution of the radiative /3~-decay mode n p+e^ +Pe+7 
with emission of a real photon |l3j-[3d| (see also |3ll. |32|). The sum of the electron-energy and angular distributions 
of these two decay modes does not suffer from infrared divergences of virtual and real photons. 

The radiative corrections to the /3~ -decay of the neutron we calculate within the standard finite-photon mass 
(FPM) regularisation of infrared divergences [l^-fl^]- As has been shown by Marciano and Sirlin the FPM 
regularisation is equivalent to the dimensional regularisation j35| . 

The radiative corrections to the lifetime of the neutron, calculated within the FPM regularisation, obey the 
Kinoshita-Lee-Nauenberg (KLN) theorem [s^. According to the KLN theorem, the radiative corrections to the 
lifetime of the neutron, integrated over the phase volume of the final state of the neutron /3~ -decay in the limit of 
the massless electron me — >■ 0, should not depend on the electron mass TOg. For the first time such an independence 
of the electron mass in the limit me has been demonstrated by Kinoshita and Sirlin [Tsj . 

We reproduce fully the radiative corrections to the lifetime of the neutron and the correlation coefficients a{Ee) 
and A{Ee), calculated in [l^-[33l a-^d [HI HI]; respectively (see Eq.Q and Appendix D). We show also that the 
correlation coefficient B{Ee) has no radiative corrections to order a/n, calculated to leading order in the large M 
expansion. This agrees well with the results, obtained by Gudkov et al. [32|. The corrections of the "weak magnetism" 
and the proton recoil we calculate in complete agreement with the results, obtained in [32| (see also [l^l). Above the 
background of these corrections the contributions of interactions beyond the SM are calculated to leading order in the 
large M expansion and in the linear approximation with respect to the Herczeg phenomenological coupling constants, 
introduced at the hadronic level in terms of the lepton-nucleon current-current interactions (see Appendix G). Such 
an approximation is good enough for the analysis of the sensitivity of the observables of the neutron /3~ -decay to 
contributions of order I0~"' beyond the SM. 

The obtained results are applied to the theoretical analysis of 1) the asymmetries Aoxp(i?e), Bc^piEe), the electron- 
proton energy distribution a{Ee,Tp), the proton-energy spectrum a(Tp) and the proton recoil asymmetry Coxp, used 
for the experimental determination of the axial coupling constant A and the correlation coefficients Aq, Bq, oq and 
Co = —xciAo + Bo), where xc = 0.27491 is a theoretical numerical factor (see section WU^ . respectively, 2) the 
lifetime of the neutron and 3) the sensitivity of the electron-proton energy distribution a(i?e, Tp), the proton-energy 
spectrum a{Tp), the asymmetries Aexp(£'e), Bcxp(Ee) and Cexp and the lifetime of the neutron Tn to contributions 
of order lO""* of interactions beyond the SM. The experimental analysis of the asymmetries Apxn(fi), BpT^ n iEp Y the 
proton-energy spectrum a{Tp) and the proton recoil asymmetry Coxp has been carried out in [lj,|37[, [381 l39l|. [40[ (see 
also [4l|) and [42|, respectively. 

II. TABLE OF CONTENTS 

The paper is organised as follows. In section Hill we calculate the electron-energy and angular distribution of the 
neutron /3~ -decay with a polarised neutron and unpolarised decay proton and electron. We take into account a 
complete set of corrections of order 1/Af, caused by the "weak magnetism" and the proton recoil, and the radiative 
corrections of order ck/tt. Then we analyse 1) in section ITVl the electron asymmetry y4cxp(i'e), which has been used 
in [l], [13 for the experimental determination of the axial coupling constant A = —1.2750(9) and the correlation 
coefficient A^^^^"^ = —0.11933(34), 2) in section IVl the antineutrino asymmetry BcxpiEe), which has been used for 

the experimental determination of the correlation coefficient Bq'^^'^^ = 0.9802(50) in [H, [s^, 3) in section IVll the 
electron-proton energy distribution a{Ee,Tp) and the proton-energy spectrum a{Tp), which may be used for the 
experimental determination of the axial coupling constant A and the correlation coefficient oq, 4) in section rVIII the 
proton recoil asymmetry Cexp, which has been used for the experimental determination of the correlation coefficient 
^^oxp) ^ -xc{Ao + Bq) = -0.2377(26) [Hi, and 5) in section [Vm] the lifetime of the neutron. In section HXl we 
propose the theoretical analysis of the sensitivity of the electron-proton energy distribution a{Ef.,Tp), the proton- 
energy spectrum a{Tp), the asymmetries Ae-^piEe)-, i?exp(£'e) and Cexp and the lifetime of the neutron to contributions 
of order 10^"* of interactions beyond the SM. We show that to linear approximation with respect to the Herczeg 
phenomenological coupling constants the Herczeg coupling constants a^^ and of the left-left and left-right lepton- 
nucleon current-current interactions (vector and axial-vector interactions beyond the SM) can be fully absorbed by 
the axial coupling constants and cannot be determined from the experimental data on the electron-proton energy 
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distribution a(Ee,Tp), the proton-energy spectrum a(Tp) and the asymmetries Acxp(£'e), ^cxp(-E'e) and CcxpiEe)- 
This agrees well with the results obtained in [4^-[4^. In the lifetime of the neutron the contributions of the Herczeg 
left -left and left -right lepton-nucleon current-current interactions (vector and axial-vector interactions beyond the 
SM) become unobservable after the renormalisation of the CKM matrix element V^d — >■ (Kid)cff = ^tid(l + ~^'^lr) 



(see section HX)) in agreement with |431-4^- In the conclusion |X] we summarised the obtained results. In Appendix 



A we give a detailed calculation of the amplitude, the electron-energy and angular distribution of the continuum- 
state /3~ -decay of the neutron to next-to-leading order in the large M expansion. In Appendix B we calculate the 
rate, the electron-energy and angular distribution of the radiative /3~ -decay of the neutron with a polarised neutron 
and unpolarised decay particles and their contributions to the rate, the electron-energy and angular distribution of 
the neutron /3~-decay. We compare our results for the branching ratios of the radiative /3~-decay of the neutron 
with the experimental data, obtained by Nico et al. [i^ and Cooper et a l. l47| ( see also [i^), and the theoretical 
results, obtained in [l^ls^]. In Appendix C within the FPM regularisation |l3|-j30l| we analyse the infrared divergent 
contributions of one-virtual photon exchanges in the continuum-state /3~-decay mode of the neutron. In Appendix D 
we give a detailed calculation of the radiative corrections, caused by one-virtual photon exchanges in the continuum- 
state /3-decay of the neutron and the radiative /3~ -decay of the neutron. We define the radiative corrections to 
the lifetime of the neutron and the correlation coefficients by the functions gn{Ee) and fn{Ee), respectively. We 
show also that in the rest frame of the neutron and in the non-relativistic approximation for the proton one-virtual 
photon exchanges induce effective scalar and tensor weak lepton-baryon couplings of order uGf/t^ and aXGp/T^, 
respectively, depending on the electron energy E^. In Appendix E we show that the radiative corrections to the 
lifetime of the neutron, defined by the function g„(i5e), satisfy the KLN theorem. This agrees with the analysis of 
the radiative corrections, carried out in [l5l |. In Appendix F we calculate the radiative corrections to the lifetime 
of the neutron, using the procedure proposed by Sirlin [l^. We confirm that the unambiguity of the observable 
radiative corrections to the lifetime of the neutron is caused by gauge invariance of the amplitude of one-virtual 
photon exchanges in the continuum-state /3~-decay of the neutron [18|. In Appendix G we calculate the contributions 
to the correlation coefficients and the lifetime of the neutron from phenomenological vector, axial-vector, scalar and 
tensor interactions beyond the SM [1]~0- Appendix H we calculate the proton recoil corrections, caused by the 
electron-proton Coulomb interaction in the final state of the neutron/3^ -decay. We show that in the experimentally 
used electron energy region 250 keV < Te ~ Ef, — me < 455 keV [1[ the contributions of these corrections to the 
correlation coefficients are of order 10^^ — 10~^ and may be neglected for the analysis of contributions of order 10^^. 
In Appendix I we give a detailed calculation of the electron-proton energy-momentum and angular distribution of the 
neutron /3~-decay with a polarised neutron and unpolarised decay proton and electron. We calculate the electron- 
proton energy distribution a(-Ee, the proton-energy spectrum a (Tp) and the correlation coefficient C of the proton 
recoil asymmetry by taking into account the 1 /M corrections from the "weak magnetism" and the proton recoil and 
the radiative corrections of order a/Tr, calculated to leading order in the large M expansion. 



III. ELECTRON-ENERGY AND ANGULAR DISTRIBUTION OF NEUTRON /3 -DECAY IN 

STANDARD MODEL 

For the analysis of the electron-energy and angular distribution of the continuum-state /3~-decay of the neutron 
we use the Hamiltonian oiV — A interactions with a real axial coupling constant A and the contribution of the "weak 
magnetism" [ssl. [sH 

Uwix) = ^ Vud[[i^p{.x)^^{l + A7')Vn(x)] + ^d^[i,j,{x)aM^)\][Mx)r{'^ - l'')^^)] (1) 

invariant under time reversal, where ipp{x), i>n{x), ipe{x) and ip^ix) are the field operators of the proton, neutron, 
electron and antincutrino, respectively, 7'', 7^ and a'"^ = ■^{'j'^^'^ — 1^1^) are the Dirac matrices [s^ and k = 
Kp — Kn = 3.7058 is the isovector anomalous magnetic moment of the nucleon, defined by the anomalous magnetic 
moments of the proton Kp = 1.7928 and the neutron k„ = —1.9130 and measured in nuclear magneton 

For numerical calculations we use Gp = 1.1664 x 10"" MeV"^ and = 0.97428(15) @. The value of the CKM 
matrix element \Vud\ = 0.97428(15) agrees well with \Vud\ — 0.97425(22), measured from the superallowed 0+ 0+ 
nuclear /^"-decays [521. I* satisfies also well the unitarity condition \Vud\'^ + IKisP + = 0.99999(43) for the 

CKM matrix elements Q. The error Ay = ±0.00043 of the unitarity condition is determined by the errors of the 
CKM matrix elements |Kd| = 0, 97428 ± 0.00015, = 0.2253 ± 0.0007 and jKbl = 0.00347t[;:[J^!^}^ (see Eq.(11.27) 
m p. 150 of Ref.H). As a result the error of the unitarity condition is equal to 



Au = |2KgAK,P = 0.00043, 
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where q = d,s,b and AKd = 0.00015, AKs = 0.0007 and AKh = 0.00016, respectively. 

The amphtude of the continuum-state /3~-decay of the neutron, calculated in the rest frame of the neutron and to 
next-to-leading order in the large M expansion taking into account the contributions of the "weak magnetism" and 
the proton recoil, is (see Appendix A) 



Gf f 

M{n pe~ De) = -2m„ --j= [we7°(l - - A [ifla fn] ■ [uel{^ ~ l^)'^v\ 

where and ipn are the Pauli spinorial wave functions of the proton and neutron and and are the Dirac bispinor 
wave functions of the electron and antincutrino, respectively. Then, is a 3-momentum of the proton related to 
3-momenta of the electron fee and antincutrino k as kp = —k^ — k, X = A(l — Eo/2M), where Eq — {rn^ — ra^ + 
TOg)/2m„ = 1.2927 McV is the end-point energy of the electron-energy spectrum, calculated for m„ = 939.5654 MeV, 
TOp = 938.2720 McV and = 0.5110 MeV Q. From Eq.© one may see that the parameter of the large M expansion 
or the 1/M corrections to the amplitude of the /3~-decay of the neutron is kp/M ^ Eq/M ^ 10~^. The detailed 
calculation of the amplitude Eq.([2]) is given in Appendix A. 

The electron-energy and angular distribution of the neutron /3^-decay takes the form P. l32l| 

- (1 + 3^^) - ^^)^ ^/W^eE.EiE., Z ^ 1) a>,, ik.rk)aE.) 

X (l + |-| + Me^) ^ + BiE,) + K^E^) • • + Qr^iE^) ikl^^ 

+ ^(i?e)^^^), (3) 

where E ~ Eq — Ee is the antincutrino energy and dfl^ and dQ are the infinitesimal elements of the solid angles of 
the electron and antincutrino 3-momenta relative to the neutron spin, respectively. 

The function (fcg, fc) is defined by the contribution of the proton recoil. To next-to-leading order in the large 
M expansion it is equal to (see Eas. (|A-19p - (|A-23p in Appendix A) 

$,-(fce,fc) = l + -^fi?e-^). (4) 



(5) 



M\ E 

The function F{Ee, Z = 1) is the relativistic Fermi function [H, [s^l (see also js^ ) 



, , 1 \ 4(2r„me/3)2^ 



2V r2(3 + 27) (1-^2)7 



r(i + 7 + «^ 



where /3 = ke/E^ = ^El — m\ / E'e is the electron velocity, 7 = — o? — 1, rp is the electric radius of the proton 
and a. = 1/137.036 is the fine-structure constant. In numerical calculations we will use = 0.841 fm [ssj . 
Following [32I we transcribe the r.h.s. of Eq.® into the form 

^^^i^;^ - + ^^'l ^ ; ^^^^ ^ - ^) 

X {1 + a(i?.) + A[E^) ^ + i?(i?.) ^ + A-„(i?e) • ^'^^l' • + Q„(i?e) ikl^^ 



FIG. 1: The densities {o./-k) gn{Ee)pp-{Ee) (left) and (a/n) fn{Ee)pp-{Ee) (riglit), measured in MeV~^, of the radiative 
corrections to the hfetime of the neutron and the correlation coefficients a{Ee) and A{Ee), where pg-{Ee) is the electron- 
energy spectrum density Ea. (|D-59|l . 



The correlation coefficients are given by (see Appendix A, B, C and D) 
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M 1 + 3A2 
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C(-Be) aiEe) = ao ( 1 + - 9niEe) + - fniEe) + 
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E, 



- 5A^ - (3k - 4) A - (k + 1) 



aEe) B{E,) = Bo 1 + - gn{E,) + — 
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M 1 + 3A2 



-2A A - (k + 1) LB, 



+ 7A^ - (3k + 8) A + (k + 1) LBe - A 



(k + 2)A+(k+1))-^ 



1 



C{E,)K„iE,)^ — 

aEe)Qn{E,) = 



1 



AI 1 + 3A2 
1 1 
U 1 + 3A2 



5A2 



{K-A)X-{K+l)^Ee, 

-{k + 2)A + (k + l)'^Eo - (TX"^ - {k + 8)X + {k + 1))E, 

aEe)D{E,)^Q, 

where the correlation coefficients oq, Aq and i?o are determined by [l| (see also (loj ' 



ao = 



1-A2 
1 + 3A2 



An = - 2 



A(l + A) 
1 + 3A2 



= -2 



A(l-A) 
1 + 3A2 ■ 



(7) 



(8) 



The radiative corrections arc determined by the functions g„(_Ee) and /„(£'e), which arc given in Ea. (jD-58|) of Appendix 
D. 

The functions (a/7r) g„(i?e) and {ol/tt) fn{Ee) describe the radiative corrections to the lifetime of the neutron and 
the correlation coefficients a{Ef.) and A{Ee), respectively. They are equal to the radiative corrections, calculated in 
[l^-llSl and [3ll. Is^. respectively. We show in Appendix D that the correlation coefficient B{Ee) has no radiative 
corrections to order a/ir. This agrees also well with the results, obtained in [s^. The densities of the radiative 
corrections (a/7r) g„(£'e)p^- (i?e) and (a/7r)/„(£'e)p^- (i?e), where p^-(£'e) is the electron-energy spectrum density 
Ea. (|D-59p . are plotted in Fig. 1. 

The coefficients Kn{Ee) and Qn{Ee) have been introduced in [Sll, and calculated within the effective quantum 
field theory, based on the hcavy-baryon Chiral Perturbation Theory (HBChPT). Our results of the calculation of the 
correlation coefficients, carried out to next-to-leading order in the large M expansion (see Appendix A), agree fully 
with the expressions, calculated in [32| . 

The correlation coefficient D{Ee) relates to a violation of time reversal invariance. The correlation coefficient D{Ee), 
calculated for weak interactions, described by the Hamiltonian Eq.([T]) with a real axial coupling constant A, is e qual 
to zero D(Ep) = 0. In the SM a non-vanishing correlation coefficient D{Ee) may appear due to long-range j56l|-|59l| 
(see also |33j ) and short-range [60j mechanisms of time reversal violation. In the long-range mechanism of time 
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reversal violation the correlation coefficient D{Ee) is induced by the electron-proton interaction in the final-state of 
the decay due to the distortion of the electron wave function in the Coulomb field of the proton [H, [H^ , the "weak 
magnetism" and the proton recoil. In the short-range mechanism of time reversal violation the correlation coefficient 
D{Ee) takes a contribution from the CP-violating phase 5 of the CKM quark mixing matrix According to [60| . 
the contribution of the long-range mechanism of time reversal violation dominates by many orders of magnitude in 
comparison with the contribution of the short-range one. As has been shown in [56| -(59l| the correlation coefficient 
D{Ee) is a function of the electron energy E^. Using the results obtained in [53, for the electron kinetic energies 
250 keV < Te < 455 keV [l, 31, Ml and the axial coupling constant A = -1.2750 we obtain that D{Ee) ^ 10"^. 
Hence the contribution of the long-range mechanism of time reversal violation to the correlation coefficient D(Ee) is 
smaller compared with contributions of order 10~^, which may be induced by interactions beyond the SM. Rece ntly 
the correlation coefficient D{Ee) has been calculated within hcavy-baryon effective field theory by Ando et al. [6ll |. 
The authors have reproduced the result, obtained by Callan and Treiman [56j . and have found a correction, which is 
smaller compared with 10^^ in the experimental region of electron kinetic energies 250 keV <Te< 455 keV hj, |37|, 1391 . 

From Eq.© the correlation coefficients under consideration, taking into account the contributions of order 1/M 
and a/TT, we define as follows 



aiE,) 
A{E,) 



Ao{l + ^UE^))+j^ 



2A(A 



(k + I)) 4A(3A 

- -C/Q 



K 



1)) 



- ao 6C{Ee), 



1 + 3A2 1 + 3A2 J 

A2-kA-(k + 1) 5A2 - (3k-4)A- (k+1) 
1 , -^0 . , „,9 E, 



1 + 3A2 



1 



1 + 3A2 



B{E,) = So , 

- BoSCiE,), 

J_ 5A2- (k-4)A- + 
M 1 + 3A2 " 

1 rA2-(K + 2)A+(K + l) ^ 

. , ■,^n Eq 



2 A (A - (k + 1)) ^ 7A2 - (3k + 8) a + (k + 
Eq H 



1 + 3A2 



-^0 SCiEe), 
1 + 3A2 E,, 



1 + 3A2 

1) ^ A2-(k + 2)A + (k+1) ml 

Ee — — 



KniEe) 
Qn{Ee) 



M 



7A2 - (k + 8) A+ (k+ 1) 
1 + 3A2 ^" 1 + 3A2 



E. 



(9) 



Using the following expansion 



^ = {l--9niEe))-SaEe), 



CiEe) 

SCiE,) 



M 1 + 3A2 
Ci = -2A(A-(«+l)), 
C2 = 10A2-4(k + 1)A + 2 
C3 = -2A(A-(ac + 1)) 



111 

(:iEo + c2E, + c3-^ 



(10) 



we transcribe the correlation coefficients a(Ee), A{Ee) and B{Ee) into the form, which is similar to that proposed by 
Wilkinson for the correlation coefficient A^^^E^) (see Eq. ipDI) '). We get 

aiE.) = .o{l + ^ (i-A2)a + 3A2) («^^° + + ) } + ^ ^"^^^O ' 



ai=4X{X^ + 1){X- {k+1)), 

02 = -26A'* + 8{k + 1) A^ - 20A2 + 8(k + 1) A - 2, 

a3 = -2A(A2-l)(A-(K + l)) 



(11) 



and 



- ^" {1 - 2A(1 + AK1 + 3A2) (^^^° + ^^^^ + ^^f ) } (1 + ? ^"(^^ 
Al = + K + {k + 2) - K X - (k + 1), 
A2 = 5X^ + K X^ - {5k + 6) X^ + 3k X + {k + I), 
A3 = -4A2(A + 1)(A- (k + 1)) 



(12) 
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and 

-""^{'-Jl 2A(l-AKl + 3,i) (^^^° + ^^^^ + }' 
Bi = -2A(A + 1)2(A-(k+1)), 

52 = A^ - (k - 4) A=^ - (5k + 2) A^ - (3k + 4) A + (k + 1), 

53 = (A2-1)(A-1)(A-(k + 1)). (13) 

For the derivation of Eq. lfTTj) . Ea. (fT2)) and Ea. (fT3)) we have neglected the terms of order {a/TT){Eo/M) ^ 10^^, which 
arc smaller compared with contributions of order 10~^ beyond the SM. 

In order to estimate the values of the obtained 1/M corrections we calculate them at A = —1.2750. This gives 

1 1 / Tfi^ \ E jn 

SCiEe) = ]g Y^^(Ci^^o + C2Ee + C3^) = -3.57 X 10-3 + 9.90 x 10"^ _£ _ 1.41 x 10"^ 

SaiEe) = jj ^3,,) {a,Eo + a.E. + ) = -3.00 x 10^ + 8.59 x 10-| + 1.41 x 10^ 

SAm . 4 ,,(,^,;,^3,.) {A^Eg + A.E. -,A<y 3.44 x 10^ + 1.47 x 10-| + 1.42 x 10- 

mE.) = -jj ^3,.) {b.E. + B.E. + B,f) = -4.66 x 10^ - 2.97 x lO-f^ + 1.36 x 10^ 

X„(£;e) = 7.14x 10-^^, 
Eo 

Qn{Ee) = 3.19 X IQ-^ - 7.27 x lO^^-^. (14) 

Due to strong dependence on the axial coupling constant A the numerical values of the contributions of the "weak 
magnetism" and the proton recoil, calculated at A = —1.2750, vary from 10"''' to 10~^ for energy independent and 
energy dependent terms. 

Summary. The correlation coefHcients of the electron-energy and angular distribution of the neutron /3~-decay with 
a polarised neutron and unpolarised decay proton and electron are calculated by taking into account a complete set 
of the 1/M corrections, caused by the "weak magnetism" and the proton recoil, and the radiative corrections of order 
a/TT, calculated to leading order in the large Al expansion. The obtained expressions for the correlation coefficients 
should be used as a theoretical background for contributions of order 10"'' of interactions beyond the SM, calculated 
in Appendix G (see Ea. (jG-8p '). Contributions of order 10~^ of interactions beyond the SM may be determined 
by measuring the asymmetries Aexp(-E'e), B^xpiEe) and Cexp between the neutron spin and the 3-momenta of the 
decay particles, the electron-proton energy distribution a{Ee,Tp) and the proton-energy spectrum a{Tp), related to 
correlations between the 3-momenta of the electron and proton, and the lifetime of the neutron t„ (see section HX)) . 



IV. STANDARD MODEL ANALYSIS OF EXPERIMENTAL DETERMINATION OF CORRELATION 

COEFFICIENT Aq. ELECTRON ASYMMETRY A^^p(Ee) 

For the experimental determination of the correlation coefficient Aq, defining correlations between the neutron spin 
and the electron 3-momentum in the SM to leading order in the large M expansion [l[ , the directions of the emission 
of the antineutrino are not fixed and one has to integrate over the antineutrino 3-momentum fc. As a result we arrive 
at the following electron-energy and angular distribution P, [s^ 

"^'^iEldnf"^ = {l + 3X')^^^{Eo-Eer,/W^eEeF{Ee,Z^l)aE,) 

X (1 + A(^)(i?e) (1 + ^ fniEe)) ' /? ) , (15) 



where we have denoted 



A{E,) + i g„(i?,) = ^(^)(i?e) (l + ^ fn{Ee)) , (16) 
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dfle = 2tt sindedOe is an infinitesimal solid angle of the electron 3-niomentum with respect to the neutron spin and 
£,n ■ (3 = P/3 cosOe with the neutron polarisation P = I'Cnl ^ 1- The correlation coefficients A{Ee) and Qn{Ee) are 
given in Eq. p^ and Eq.®, respectively. The contribution, proportional to Qn{Ee), with structure ^„ • ke/Ee we 
obtain having integrated the term with the structure ■k){ke-k) / E^E"^ in Eq. ^ over directions of the antineutrino 
3-momcntum k. 

The asymmetry, which may be used for the experimental determination of the axial coupling constant A and the 
correlation coefficient Aq^ takes the form 

where N^{Ee) are the numbers of events of the emission of the electron forward (+) and backward (— ) with respect 
to the neutron spin into the solid angle Ar2i2 ~ 2tt{cos9i — cos 6*2) with < (p < 2tt and 9i < dg < 62- They are 
determined by [62|] 

N+{Ee) = 2^N{Ee) ' (i + A''^\E,) (1 + ^ fn{E,S) PP cos^,,) sin^ed^e = 
= 2t:N{E,)(i + i {E,) (l + ^ /„(Se)) Pp (cos 0i + cos 62)) (cos - cos 62), 

/7T — O2 
{1 + A'^^HE,) (1 + ^ .fn{Ee)) PP COS^e) Siu d^e = 

= 2t:N{E,)(i ~ i A(^) (£;e) (1 + ^ /n(Se)) PP (cos 0i + cos 62)) (cos 01 - COS 62), (18) 
where N{Ee) is the normalisation factor equal to 

N{E,) = (1 + 3A2) - E,f ^El - ml E, F{E,,Z = 1) ({Ee). (19) 

The correlation coefficient A'-^'>{Ee) is 

4"^^ = ^ ( - SA-'^ + (Sk + 5) A2 - (2k + 1) a - (k + 1)) = -2(a - {k + 1)) (^A^ - ^ A - i) , 



4"^ 



( - SA"^ + {3k + 12) A-'' - (9k + 14) A^ + (5k + 4) A + (k + 1)) = -2(a - (k + l)j (a^ " + ^ ^ + ^) ' 
4A2(A + 1) (a- (k + 1)). (20) 



It agrees well with the result, obtained by Wilkinson [sj, |6^ . We note that the correlation coefficient A{Ee) + ^Qn {Ee) 
differs from the Wilkinson correlation coefficient A'^'^\Ee) by the contribution of the radiative corrections, described 
by the function {a/i:) fr,{Ee). In the replacement ^(£'e) + iQ„(i;e) A'^'^\Ee){l + {a/T:)fn{Ee)) we have neglected 
the contributions of order {a/'K){EQ/M) ^ 10^^, which are smaller compared with contributions of order 10^"* of our 
interest. The contribution to the correlation coefficient A'^'^^E^) of the "weak magnetism" and the proton recoil, 
calculated at A = —1.2750, is equal to 

= -8.55 X 10"^ + 3.73 x 10"^ ^ ^ 42 x IQ-^ (21) 

Eo Ee 

In the experimentally used region of electron kinetic energies 250 keV < Tg < 455 keV [l], [13, HI] the radiative 
corrections (a/7r) /„(i?e) vary over the region 1.53 x 10"'^ > {a/ir) fn{Ee) > 1.04 x 10^"^ and increase the absolute 
value of the correlation coefficient A^^\Ee). 

Summary. The Wilkinson expression for the correlation coefficient A^'^^E^) (sec Eq. (|^D|) ). taking into account a 
complete set of the 1 /M corrections from the "weak magnetism" and the proton recoil, is improved by the account for 
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the radiative corrections, described by the function 2.81 x 10"'^ > (a/7r) fn{Ee) > 0.62 x 10~^ for the electron energies 
nie < Ee < Eq (sec Ea. (jD-58|) ). The contribution of the radiative corrections increases the absolute value of the 
correlation coefficient The results, obtained in this section, should be used as a theoretical background for 

the experimental determination of contributions of order 10~^ of interactions beyond the SM from the experimental 
data on the electron asymmetry A^xpiEe) (see section HX]). 



V. STANDARD MODEL ANALYSIS OF EXPERIMENTAL DETERMINATION OF CORRELATION 
COEFFICIENT Bq. ANTINEUTRINO ASYMMETRY Bexp(£e) 

In the SM to leading order in the large M expansion the correlations between the neutron spin and the antineutrino 
3-momentum are defined by the correlation coefficient Bq [i|, which may be determined from the experimental data 
on the asymmetry Bcxp{Ee), defined by [s^ 

N-iE,)^N++iE,) 
^cM^e) - N—{E,) + N++iE,)- ^ ' 

It defines the asymmetry of the emission of the antineutrinos into the forward and backward hemisphere with respect 
to the neutron spin, where N^^{Ee) is the number of events of the emission of the electron-proton pairs as functions 

of the electron energy E^- The signs (++) and ( ) show that the electron-proton pairs were emitted in (++) and 

against ( ) a direction of the neutron spin. This means that antineutrinos were emitted against (++) and in ( ) 

a direction of the neutron spin. The number of events N (Ee) and N^^{Ee) are defined by the electron-energy 
and angular distribution of the neutron /3~-dccay, integrated over the forward and backward hemisphere relative to 
the neutron spin, respectively. 

The integration region for the electron-proton pairs, emitted in the direction (H — h) of the neutron spin, is defined 



by the constraints [6j|: ^„ • fcg = Pke cos6e > and ^„ ■ kp = ^„ • (—fee ^ k) = PE{—r cosOe — cosfl) > or 



—rcos^e > cos 6', where r = k^/E = ^ El — ml/ {Eq — Ee) and P = |^„| < 1 is the neutron polarisation. For 
N^'^{Ee) and r < 1 we obtain the following expression 

N++{Ee) = 27riV(ii;e){ (l -\ap + _ i A'„ /?) - ^p(^B - ^ Q„/3 

- (l + |PA/?) +K'(al3 + ^PB + ^P /i„ (5^^ - l^r^ P Q„ i3 - lao r{l - r^) ^} , (23) 

where N{Ee) is the normalisation factor Ea. (jl9p . For r > 1 the upper limit of the integration over cos^e is restricted 
by cos6'e < l/r. The result of the interaction is 

11/. 2„„\ 11/ „ 4„.„ 4„„ „\ 1 1 1 ^-^If. 

r\ r^J Ml' 
(24) 



N++{Ee) = 2^N{Ee){--(l ^Pb) --(ap - -PA(3 - -PQ„/3) - j^^PKnfS' + -aol3' 



For the calculation of N (Ee) we have to integrate over the region [6J]: ■ fee = Pke cos0e < and £,„ ■ kp — 
PE{—r cosde — 008 6*) < or cos9 > —rcosOe- The number of events N {Ee) for r < 1 is given by 

N—{Ee) = 2T,N{Ee)[(l-\aP-]^Pfi[A-l^K,,p) + Iq,,/?)) 

-^r(l-^PA/3) +ir2 (a/3-1 PS- lp/C„/32) +lr3pQ„/3-iao/3'r(l-r2)|^}. (25) 

For r > 1 the lower limit of the integration over cos 6*2 is restricted by cos 6*6 > — l/r. The number of events N (-^e), 
calculated for r > 1, is equal to 

N-[Ee) = 2.N[Ee){\\{l + Ipb) - ll(a/3 + \pAp + \pQ..p) + + ^" 7 - Si' 

(26) 

Using our formulas for the numbers of events we calculate the asymmetry B^xpiEe)- For r < 1 or < 7e < 
{Ea - me)'^/2Eo = 236 kcV and for r > 1 or (Eq ~ me)^/2Ea = 236 keV < Te < Eq ~ nie the asymmetry B^xpiEe) is 
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T, [keV] 



FIG. 2: The antineutrino asymmetry -Bcxp(-Ee), including a complete set of the 1/M corrections, caused by the "weak mag- 
netism" and the proton recoil, and radiative corrections of order q/tt, calculated to leading order in the large M expansion; 
Te — — rrie is the electron kinetic energy. 



equal to 



i3exp(i?e) = ^ (S(3 - r2) - (3 - 2r)Af3 + (l - ^r^) K^13^ - (l - ^r^^Q^/S ) 
((8 - 4r) - (2 - r')a(5 - a^P^ r(l - r^) ^) '\ 



and 



(27) 



(28) 



respectively. At r = 1 or Tg = (i?o — m^Y j^E^ ~ 236 kcV the asymmetry i?exp(^'e) is continuous. To leading order 
in the large M expansion the asymmetry Ba-xp{Ee) reduces to the form 



2P 



where the correlation coefficients B, A and a are equal to 

Bo 





- 


(1- 


)A/3 


(2-.)4 


('- 


-'r1 


al3 



r <l 



r > 1, 



B 



A 



Aa[l + -fn{Ee 

n 



B-'-A^- 
2 r 

1 - - a - 

4 r 



ao l + -/„(£^e 

TT 



(29) 



(30) 



In Fig. 2 we plot the asymmetry Bo^piEe), given by Ea. (P7)) and Eq. (P5|) and obtained in the SM with the account 
for the contributions of the 1/M corrections, caused by the "weak magnetism" and the proton recoil and the radia- 
tive corrections of order a/n, calculated to leading order in the large M expansion and described by the function 
(ct/n) fn{Ee). In [ssl. [s^ for the experimental determination of the correlation coefficient Bq the experimental data 
on the asymmetry B^xpiEe) were fitted by Ea. (|29p with the replacement B ^ Bq, A ^ Aq and a — > ao, respectively. 
Such an asymmetry was calculated in [6J] . 

The theoretical value Bq = 0.9871(1) of the correlation coefficient Bq, calculated for A = —1.2750(9), agrees 1) 

within 1.5 standard deviations with the experimental values S^^''^^ = 0.9802(50), B^^''^' = 0.9801(46) and B^"""^'^ = 
0.9894(83), obtained in [H, [s^ (see also [l|), [i^ and [6^, respectively, and within two standard deviations with the 
experimental one b('=''p) = 0.967(12), obtained in^. _ „ „ 

Summary. The antineutrino asymmetry Bg^piEe) 

(oxp) ^ 



determination of the correlation coefficient Bn 



calculated in [6j| and used in [38|, |39[ for the experimental 
0.9802(50), is improved by the contributions of a complete 



set of the 1/M corrections, caused by the "weak magnetism" and the proton recoil, and the radiative corrections 
{a/ir) fn{Ee), calculated to leading order in the large M expansion. The obtained expression of the antineutrino 
asymmetry BcxpiEe) should be used as a theoretical background for the experimental determination of contributions 
of order 10^^ of interactions beyond the SM (see section HX)) . 
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VI. STANDARD MODEL ANALYSIS OF EXPERIMENTAL DETERMINATION OF CORRELATION 
COEFFICIENT oq. ELECTRON-PROTON ENERGY DISTRIBUTION a{Ee,Tp) AND PROTON-ENERGY 

SPECTRUM a{Tp) 

For the experimental determination of the correlation coefficient oq the correlation coefficient a{Ee), given by 
Eg. dllD . can be hardly used, since the antineutrino is hard to detect. Thus, for the determination of oq and the 
contributions of interactions beyond the SM one should measure correlations of the 3-momenta of the decay charged 
particles. Using the results, obtained in Appendix I, the electron-proton energy spectrum for the neutron /3~-decay 
with unpolarised particles can be given in the following form 

^^dEAT^"^ = (t„)smM (1 + 3A2)^M^ a{E,, T^) (l + ^ g^{E,)) F{E,, Z^l) E,, (31) 

where (r„)sM is the theoretical lifetime of the neutron, calculated in the SM (see section rVIII[) and Tp is the kinetic 
energy of the proton varying from zero to its maximal value (Tp)niax = {f^n — rup)'^ — mg/2m„ = {Eq — m'l)/2M = 
0.751 keV, i.e. < Tp < 0.751 keV. The limits of the integration over the electron energy Ee, i.e. (i?e)min < ^-e < 
(i?e)max, are the functions of the proton kinetic energy Tp. They are adduced in Appendix I. The electron-proton 
energy distribution a(Ee,Tp) is defined by (see Appendix I) 

aiE,,Tp) = CiiE,,Tp)+ao (l + + f f^i^-)) UEe,Tp). (32) 

The functions C^i{E(,,Tp) and C,2{EeiTp) are given by in Appendix I. They are calculated to leading order in the 
large M and do not contain the radiative corrections. In the electron energy region (£'e)min < < (i?e)max the 
contribution of the radiative corrections {ol/tt) fn{Ee) to the electron-proton energy distribution a{Ee,Tp) relative to 
the correlation coefficient ag is of order 10~^. The account for the radiative corrections is very important for a correct 
experimental determination of contributions of order 10""* of interactions beyond the SM. 

Integrating the electron-proton energy spectrum Eg. pip over the electron energy (i?e)min < E^, < (£'e)max we 
obtain the proton-energy spectrum 
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where a{Tp) is defined by 



^j.^ -(r„)sMM(l + 3A^) -^;;;3^l a{Tp), (33) 



a{Tp) - .91 (Tp) + ao (l + ^) 92{Tp). (34) 



The functions gi(Tp) and 52 (Tp) are given in Appendix I. Recently the correlation coefficient ag and the axial coupling 

constant A have been determined from the proton-energy spectrum by Byrne et al. [10] . The obtained value a'^^^'^ = 

-0.1054(55) can be fitte d by the axial coupling constant A — —1.271(18). In turn, the experimental value Oq — 
—0.1017(51), obtained in [41|, defines the axial coupling constant equal to A = —1.259(17). These experimental values 
agree with the theoretical value of the correlation coefficient ao = —0.1065(3), calculated at A = —1.2750(9), and 
with the axial coupling constant A = —1.2750(9) within one standard deviation. 

From the point of view of the experimental determination of contributions of order 10~^ of interactions beyond the 
SM (see a discussion in section VII) the experimental accuracy of the determination of the axial coupling constant 
A by measuring the proton-energy spectrum as well as the electron-proton energy distribution a{Ef,,Tp) should be 
improved by more than two orders of magnitude in comparison with the accuracy of Byrne's experiment [40j . There 
are two major, funded experiments, which are currently attempting to do this. They are 1) aCORN at the National 
Institute of Standards and Technology (NIST) and 2) Nab at the new Spallation Neutron Source (SNS) in Oak Ridge 
at Tennessee [68| . The expected experimental accuracy of the determination of the correlation coefficient ag in the 
aCORN and Nab experiments is better than 1 % (see also [6^ ) . 

We have discussed the experimental determination of the correlation coefficient oq by measuring the electron- 
proton energy distribution a(Ee,Tp) and the proton-energy spectrum a(Tp). However, there are another possibilities 
to determine experimentally the correlation coefficient gq. It can be extracted from the experimental data on the 
(TejCOS^ep) and (Te,cos0ep) distributions [70| . where Oen and 6ep are angles of the electron-antineutrino and electron- 
proton correlations, respectively. 

Summary. The electron-proton energy distribution a{Ee,Tp) is calculated by taking into account a complete set 
of the 1/M corrections, caused by the "weak magnetism" and the proton recoil, and the radiative corrections of 
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order a/7r, calculated to leading order in the large M expansion. The proton-energy spectrum is improved in 
comparison with that, used in [40j, by the account for the 1/M and radiative corrections. The obtained expression 
for the electron-proton energy distribution a{Ef,,Tp) as well as the proton-energy spectrum a{Tp) should provide 
a theoretical background for the experimental determination of contributions of order 10~* of interactions beyond 
the SM (see section HX)) . For the analysis of contributions of interactions beyond the SM we propose to multiply the 
electron-proton energy spectrum of the rate of the neutron /3~ -decay by the lifetime of the neutron t„, calculated 
with the account for the contributions of interactions beyond the SM (see Appendix G). As a result the analysis 
of the contributions of vector and axial-vector interactions beyond the SM by means of the electron-proton energy 
distribution r„a(£'e, Tp) and the proton-energy spectrum Tna{Tp) reduces to the analysis of these contributions to the 
axial coupling constant A only (see section HX)) . 



VII. STANDARD MODEL ANALYSIS OF EXPERIMENTAL DETERMINATION OF CORRELATION 

COEFFICIENT Co. PROTON ASYMMETRY Ccxp 

The correlations between the neutron spin and the proton 3-momentum are described by the correlation coefficient 
C (zli-izEl- It defines the proton recoil asymmetry Co — —xc{Ao + Bq), where xc is the theoretical numerical 
factor, calculated within the SM to leading order in the large M expansion. The first measurement of the correlation 
coefficient Cq™^'' = —0.2377(26) has been performed by Schumann et al. The angular distribution of the 

probability of the neutron /3~-decay, related to the proton recoil asymmetry, is given by [7l| 

4^^^I^lM = i + 2PCcos0p, (35) 

Cii Lip 

where (Kip = 2tt sm9pd9p is the infinitesimal solid angle of the proton 3-momentum with respect to the neutron 
polarisation i.e. ■ fcp = PkpCosOp and P = |^„| is the neutron polarisation. The correlation coefficient C, 
calculated at the account for the 1/M and radiative of order a/n corrections, is equal to (see Appendix I Ea. (|I-29p ) 

„ 1 Xr , , „ , 1 Xq , a 1 X-[ n , a 1 X-\ i „ 1 1 / , 1 X-\ 2 , , , 1 X^ -i 

^ = -2X7(^° + ^")+2]^^" + ^2if^°-^2i'^°+MTT3A^(^2i'-('^ + l)^2^ 

1X14 , ^^ 1^15 + i"3 , .,lXi6 + l4\ alX2 1 1 f^X. ^ 



" J^TT^ - " - 5^Tr^ + + + \y[ - - + 1'"' if)- 

where the numerical factors Xj (j = 1, . . . , 14) and Yj (j = 1, 2, 3, 4) are calculated in Appendix I (see Ea. (|I-22p and 
Eg. (1^26)) ). The factor xc = Xs/2Xi = 0.27491 agrees well with the factor xc = 0.27594, calculated by Gliick 
As a result the correlation coefficient Co is equal to 

Co = -0.27491 (^0 + Bo). (37) 

The numerical value Co = —0.2386, calculated at A = —1.2750, agrees well with the experimental one Cq^^ — 
—0.2377(26) [1^. Defining the proton recoil asymmetry Coxp as we have defined the asymmetry Acxp{Ee) (see 
Eq.dlTl)) we obtain 

Coxp = PC (cos?9i + cosz?2), (38) 

where the polar angles 9i and 02 define the solid angle Aili2 = 27r(cos6'i — 00362) of the proton emission to the 
forward and backward hemisphere with respect to the neutron spin. 

Summary. The correlation coefficient C, describing correlations between the neutron spin and the proton 3- 
momentum, is calculated by taking into account a complete set of the 1/M corrections, caused by the "weak mag- 
netism" and the proton recoil, and the radiative corrections of order a/n, calculated to leading order in the large 
M expansion. The obtained result should be used as a theoretical background for the experimental determination of 
contributions of order 10~* of interactions beyond the SM (see section HXl) . 



VIII. STANDARD MODEL ANALYSIS OF LIFETIME OF NEUTRON 

Having integrated the electron-energy and angular distribution Eq.® over the directions of the electron 3- 
momentum fcg, the antineutrino 3-momentum k and the electron energy within the limits me < E^ < Eq, 
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we obtain the rate of the neutron /3 -decay. It is equal to [T3| 



(AJsM = (l + 3A2)^M^/„(i?o,Z = l), (39) 



where the Fermi integral fniEa, Z = 1), given by 



MEo,Z^l)^ / {Eo - E,f ^El ~ ml E, F{E,, Z = 1) 1 + - gr,{E,) 



X {l + ^ [(^"^' - 4(«: + 1) A + 2) E, -2X(^X-iK + 1)) (Sq + ^)] } dii;e, (40) 

contains the contributions of the "weak magnetism" , the proton recoil and the radiative corrections, described by the 
function gn{Ee). The calculation of the lifetime of the neutron with the Fermi function, determined by Eq.(IS|), the 
axial coupling constant A = —1.2750(9) and the CKM matrix element Vud = 0.97428(15), gives (T„)gM = 879.6(1.1) s. 
The error bars ±1.1 s are defined by the error bars of the experimental value of the coupling constant and the CKM 
matrix element. The theoretical value of the lifetime of the neutron (t„)sm = 879.6(1.1) s agrees well with the 
experimental values Tn^^'' = 878.5(8) s and Tn^^'' = 880.7(1.8) s and Tn^'^'^ = 881.6(2.1) s, measured by Screbrov et 
al. [l2| . Pilchmaier et al. [76l | and Arzumanov et al. fll\ . respectively. It agrees also well with the new world average 
values (w.a.v.) of the neutron lifetime r,!'^ ''"-^ = 880.1(1.1) s Q, r^'^ ''"-^ = 880.0(9) s ^ and r^'^ ''"-^ = 881.9(1.3) s 
(79j . respectively. 

Summary. The lifetime of the neutron is calculated by taking into account the radiative corrections by Sirlin et 
ai, calculated to leading order in the large M expansion, and a complete set of the 1/M corrections, caused by 
the "weak magnetism" and the proton recoil. The obtained result should be used as a theoretical background for 
the experimental determination of contributions of order 10~^ of interactions beyond the SM (see section HXj) . The 
numerical value of the lifetime of the neutron (t„)sm = 879.6(1.1) s, calculated in this section, is left unchanged after 
the absorption of the Hcrczcg phenomenological coupling constants a^^ and a^^ of the left-left and left-right lepton- 
nucleon current-current interactions (vector and axial-vector interactions beyond the SM) by the axial coupling 
constant A — ?> Xcs = (A — + a^^) /(I + a^^ + a'lj^) = A — a^^ + a'lj^ — X{a\^ + a^^^) and the CKM matrix element 
Vud (Kd)off' = Kd(l + aJ^L + a^LR) (see section HX]) . 



IX. SENSITIVITY OF ELECTRON-PROTON ENERGY DISTRIBUTION, ASYMMETRIES Aoxp(Se), 
B^^p{Ee) AND Coxp AND LIFETIME OF NEUTRON r„ TO CONTRIBUTIONS OF ORDER 10^ OF 

INTERACTIONS BEYOND STANDARD MODEL 



In this section we propose a theoretical analysis of the sensitivity of the electron-proton energy distribution 
a{Ee,Tp), the proton-energy spectrum a{Tp), the asymmetries Aexp(i'e), B^xpiEe) and Coxp and the lifetime of 
the neutron t„ to contributions of order 10""' of interactions beyond the SM. 

The electron-proton energy distribution a{Ee, Tp) and the correlation coefficients a{Ee), A{Ee) and B{Ee), including 
the contributions of interactions beyond the SM, can be written in the form 

a{Ee,Tp) = Ci{Ee,Tp),tf + a,s{E,)(l + ^f„{Ee))C2{Ee,Tp), 

a{Ee) = a,s{Ee)(l + ^fn{Ee)), 
A{E,) = Aeff(-Be)(l + ^/n(£^e)), 

B{Ee) = Seff(^e), (41) 

where Ci(^e, Tp)cff , (icsiEe), aoff(^'e)j ^off(^'e) and Bcs{Ee) are the electron-proton energy distribution and the 
correlation coefficients, defined by the contributions of the "weak magnetism" , the proton recoil and interactions 
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beyond the SM only. They are given by 

11 / T1~L \ 

1 1 / nig \ nie 

ac«(Ee) = (ao)ofr + — , ^^n,^y [aiEo + a2Ee + j + "^4 — , 

Aeff(Se) = [AoU + ir J^^TTT^::^ (AlEo + A2E, + As^j +^4"^, 

2 

Seft(£;e) = iBo)eS + ij T^^T^ ( ^iEq + B^E, + B^^^ + Bi^, (42) 

where {rrie/ Ee)sM is the average value, calculated in the SM with the electron-energy spectrum density Eq. (jD-59[) . 
The terms, proportional to {rrie/ Ee)sM, come from the lifetime of the neutron, calculated at the account for the 
contributions of interactions beyond the SM (see Appendix G). The coefficients aj, Aj and Bj for j ~ 1,2,3 are 
given in Ea. (fTT|) . Ea. ([T2|) and Ea. (fT3)) . respectively. They are calculated to next-to-leading order in the large M 
expansion and include the contributions of the "weak magnetism" and the proton recoil only. The numerical values 
of these corrections are estimated in Eq. p^ at A = —1.2750. The terms {ao)cS, (Ao)off and (Bo)ofr arc the sums of 
the correlation coefficients ag, Aq and Bq and energy independent contributions of interactions beyond the SM (see 
Appendix G). They read 

(ao)off = ao+ (^^^\^2y Rc{SCv - SCy) + 4XRc{5Ca - 6Ca) 

{AoU = Ao+ ^^^3^,^, ( - A(3A2 - 2A - 1) Re{SCv - SCy) - {3X^ - 2A - l)Re{6CA - SCa) 
(Bo)eff = Bo+ ^^^\yy_ ( - A(3A' + 2A - 1) ReidCv - SCy) - (SX'' + 2X-1) Re{5CA - SCa) ) • (43) 
The coefficients 04, A4 and B4 are induced by interactions beyond the SM only. They are equal to 

«4 = -7T—:^(MCs - Cs) + 3ARc(Ct - Ct) 



Aa^2 



(1+3A2)2 

A(l + A) 



(Re(Cs - Cs) + 3ARe(CT - Ct) 



(1 + 3A2)2 

B^ = ^Yl + 3A2^)2^^^ (^^'^^^ " ~ ^'^^^ " ■ ^^^^ 

For the derivation of Eqs.(H5]) and we have used Ea. (|G-8|) in Appendix G with the coupling constants Cy = 
1 + SCy, Cy = — 1 + SCy, Ca = —A + SCa and Ca = A + SCa and have kept only the linear contributions of the 
deviations from the coupling constants of the SM. By assumption the axial coupling constant A in Eas. p5)) and (|31]) 
is determined by all interactions within the SM only. 

The differences of the phenomenological coupling constants SCy — SCy, SCa — SCa, Cs — Cs and Ct — Ct, calculated 
in terms of the Herczeg phenomenological Icpton-nucleon coupling constants 0] (see Appendix G), take the form 

5Cy-SCy = 2{alL + aU), 
SCa -SCa = 2{a1^-a1j,), 

Cs-Cs = 2{Al^ + Alj,) = 4as, 

Ct-Ct = 4:a'l^^4:aT- (45) 

This means that the phenomenological coupling constants SCy — SCy, SCa — SCa, Cs — Cs and Ct — Ct describe 
interactions beyond the SM of left-handed leptonic currents and left (right) -handed hadronic currents, i.e. L(E}L and 
L (E) R, respectively. 

The experimental determination of the correlation coefficients ag, Aq and Bq runs as follows. Using the theoretical 
expressions for the electron-proton energy distribution a{Ee,Tp) and the asymmetries Ac-xp{Ee) and B^xpiEe) the 
experimental data on them are being fitted by a tuning of the axial coupling constant [l| . Such a procedure gives the 
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axial-coupling constants Aq, Xa and \b, obtained from the fit of the experimental data on a{Ee,Tp), Acxp{Ee) and 
BcxpiEe), respectively. In terms of these axial coupling constants one may define the correlation coefficients (ao)^^^\ 
(Aq)^^^^ and {Bo)l^g'^\ respectively, as follows 

f s(exp) _ 1 - ...(cxp)_ Aa(1 + Aa) .p,^(oxp)_ Ab(1 - \b) , . 

The theoretical expressions for the axial coupling constants Aq, Xa and Xb in terms of the axial coupling constant A, 
defined by interactions within the SM only, and the contributions of interactions beyond the SM are 

Xa = X-^(^XRe{SCv -6Cv)+Re{SCA-SCA)y 
Xa = X-^(^XRe{SCv -SCv)+Re{SCA-SCA)), 

Xb = A - i (^XRe{SCv - 6Cv) + Re{SCA - SCa)) ■ (47) 

Thus, in the linear approximation with respect to the deviations of the phenomenological coupling constants of 
interactions beyond the SM from the coupling constants of the SM we obtain that Aa = Aa = As. This agrees well 
with the results obtained in [islj-psj. Replacing Aa, A^ and Xb by AeS, which includes the contributions of vector and 
axial-vector interactions beyond the SM in addition to the contributions of interactions within the SM, and denoting 

bp = , , (MCs -Cs)+ 3AeffRe(CT - Ct)) - ^ , l,^ f Re(as) + SAch Rc(aT)) (48) 



and 



"^ST = 1 , (AcffRc(Cs - Cs) - Re(CT - Ct)) = —^-^(x,sRc{as) - Rc(aT)), (49) 



where hp is the Fierz term (see Appendix G), we obtain the electron-proton energy distribution a(Ee,Tp) and the 
correlation coefficients a{Ee), A{Ee) and B{Ee) in the following form 

a(i^e, T,) ^ (l - bp{^)J { (l + 6.^) C^[E.,Tp) + (l + ^ + ^ /"(^^)) C.(i^e, T,)}, 

am ^ «o(l-5.^) (l + ; UE.^)[l + ^,J^ ^ 3^,^^ {am + a,E. + a,^) }, 

= H'-'^^-i) (l + ; {l - ]^ 2Ae.(l + a2)(1 + 3A.^,) (-^^^° + + -^-f ) }' 

^ - ^Atf/- A:t^^ --i) - ^7 .A.(l - A:)(1 . 3A,.) (^--0 . B.E. . E.f)}, (30) 

where the coefficients Uj, Aj and Bj for j = 1, 2, 3 are defined in Ea. (fTT|) . Ea. ([T^ and Ea. (fT5|) . respectively, with the 
replacement A — > Agft- The same replacement defines the correlation coefficients ao, and Bq in terms of Aeff (see 
Eq.®). 

The correlation coefficients a{E(.). A{Ec.) and B{Ee), extended by the contributions of interactions beyond the SM, 
together with the correlation coefficients Kn{Ee) and Qn{Ee), caused by the contributions of the 1/M corrections 
from the "weak magnetism" and the proton recoil only, determine the asymmetries ^exp(£'e) and i3exp(£'e)- For the 
calculation of the electron asymmetry A^xpiEe), taking into account the contributions of interactions beyond the SM, 
we have to replace the correlation coefficient A'^^^Ee) in Eg. pT)) by the expression 

^ ^ V "^EJX a/ 2Acff(l+Aeff)(l + 3A2jj) V 1 2 c 3 ^J/- V ) 

where the coefficients are given in Eq. ([201 with the replacement A Acff . 

The antineutrino asymmetry J5exp(^'e) is defined by Ea. (P7|) and Ea. (P5| for r < I and r > 1, respectively. For 
the account for the contributions of interactions beyond the SM the correlation coefficients a{Ee), A{Ee) and B{Ee) 
should be taken in the form, given by Ea. (j50p . 
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The proton recoil asymmetry Cexp completes the set of asymmetries, which can be measured in the neutron /3~- 
decay with a polarised neutron and unpolarised proton and electron. The correlation coefficient CeS (see Appendix 
I), defining the asymmetry Coxp and extended by the contributions of interactions beyond the SM, takes the form 



n (a o a , A(l + Acff ) , Xi7 1 (l + Aeff)(l~3Aeff) ^18 



where Csm = C + xc{Ao + Bq) and xc = 27591 (see Eq.((57)) '). The numerical factors Xir/Xi = -0.90187 and 
Xis/Xi = 0.39806 are calculated in Appendix I (see Eq. p^ ). 

Using the results, obtained in Appendix G, the theoretical expression for the rate of the neutron /3~-decay, including 
the contributions of interactions beyond the SM, is 

A„ = {Xn)sM(l + Y^^[ReidCv~SCv) + i\Re{6CA-SCA))+bF{'^)^J = 
27r3 



UEo.Z = 1) (1 + 3A2) (l + {Re{5Cv - SCy) - 3ARc(5Ca - 5Ca)) + ^J^(^)g„) - (53) 



where (A„)sm is the lifetime of the neutron, calculated within the SM (see Eg. ([5^ and Ea. ipO)) ). and {me/Ee)sM is 
the average value, calculated with the electron-energy spectrum density Ea. ()D-59p . Now we have to define the rate 
of the neutron /3~-decay Eq. (j53p in terms of the axial coupling constant Aoff , which is related to the axial coupling 
constant A as 

A = Aeff + i(AeffRc(,5Cv - SCv) + Rc{SCa - SCa)). (54) 

Substituting Ea. (|M)) into Eg . ([55)) and keeping only the linear terms in power of Rc{SCv — SCv) and Rc{SCa — SCa) 
one may show that 

(1 + 3A2) (i + -^-^(Rc(,5Cy-JCy)-3ARe(JCA-<5CA)) +foF(^)^J = 

1 + I Rc{SCv - SCv) (1 + 3A,\) (l + bF{^^)^J ■ (55) 

This gives the rate of the neutron /3~ -decay equal to 

= ^^j^f'' fn{Eo, Z = 1) (l + i Rc{SCv - SCv)) (1 + 3A,\) (l + bF{'-^)^J , (56) 

where the Fermi integral is given by Eg. ipp]) with the replacement A — )- Aoff. 

If we introduce again (A„)sm, defined by Eo. ((39| and Eg. (|40)) with the replacement A — > Acs, we get the following 
expression for the rate of the neutron /3~-decay, corrected by the contributions of interactions beyond the SM taken 
to linear approximation with respect to the Herczeg phenomenological coupling constants 

A„ = (A„)sm (l + Reial^ + a'l^)) (l + bF{^)^J , (57) 

where we have set Re{SCv — SCv) — 2Re(a^^ + ^lr)- Thus, at first glimpse a deviation of the experimental values 
r„ = 1/A„ of the lifetime of the neutron considered relative to the theoretical value of the lifetime of the neutron 
(tvOsm = l/(A„)sM, calculated in the SM at zero Herczeg coupling constants, may give an information about the 
contribution of the Herczeg left-left and left-right lepton-nucleon current-current interactions (vector and axial-vector 
interactions beyond the SM) by using the experimental value of the Fierz term bp, determined from the experimental 
data on the electron-proton energy distribution a{Ee,Tp), the proton-energy spectrum a{Tp) and the asymmetries 
^oxp(£'e), BcxpiEe) and Cexp- However, the problem is in the following that the Herczeg phenomenological interactions 
beyond the SM, when they exist, should exist always together with the interactions of the SM, and the decomposition 
of these interactions is rather artificial and based on the smallness ones with respect to others. 

Hence, using the Hamilton Eg. (|G-ip and Ea. (jG-2p we may redefine the axial coupling constant and the CKM 
matrix element as follows AoS = (A — a^^ + a'lj^)/{l + a'lj^ + a'lj^) and (Kid)off = Vud{^ + o-'ll + ^la) i* l^as been 
proposed in [4^-[4^. After such a change one may show that the rate of the neutron /3~-decay may contain only the 
contribution of the Fierz term 



An 



(A„)sM(l + fo^^(^) ), (58) 

V \ tip I SM/ 



18 



where (A„)sm is given by Eg. ip^ and Ea. (|i(7)) with the replacement A Aoff and Vud — >■ {Vud)eS- 

The definition of the efi'ective couphng constant Aeff = (A — a^^ + a'lff)/{l + a\j^ + aj;^) and the CKM matrix 
element (Kid)eff = Kid(l + ^ll + '^la) the Hamiltonian level introduces the imaginary parts to the axial coupling 
constant and the CKM matrix element 

Aoff = Re Aoff + ilm Aoff, 

(Kd)cff = Kd(l+Rc(aL+«D)(l+*M4L+«D) (59) 

where ImAoff is equal to 

ImAoff = -(1 + ReAoff)Im(a^^i) + (1 - ReAoff)Im(a^fl). (60) 

An information about the imaginary part of the axial coupling constant Aoff one may obtain by measuring the 
correlation coefficient D{Ee), describing a violation of time reversal invariance. From Ea. (|G-8[) and Ea. (|M)) we obtain 

D{E,) = DsuiEe) - Y^|^(AoffIm(aii + alj,) + Im(aii - a'lj,)) - 

= Z?SM(i;e)- Y^|^((l + Aoff)Im(a^J-(l-Aoff)Im(a^^)) -i?SM(i?e) + |^^^, (61) 

where we have replaced ReAoff by Acg, having neglected the contribution of the imaginary part, and DsM{Ee) is 
the contribution to the correlation coefficient D{Ee), calculated within the SM [5^-[6l|. As we have estimated in 
section Hill such a contribution, caused by the electron-proton interaction in the final state [s^-js^ [6l[, is of order 
10^^ for the electron kinetic energies 250 keV < Te < 455 keV. Thus, the correlation coefficient D{Ee), defined in the 
same energy region, should be sensitive to the contributions beyond the SM of 10~^. Recently the experimental value 
Dcxp{Ee) = (— 4 ± 6) X 10~^ of the correlation coefficient D{Ee) [HQ has been substantially improved with a result 
Dcxp{Ee) — (— 0.96± 1.89stat ± 1-Olsyst) X 10~^ [13 ■ However, the new experimental value as well as the old one still 
implies that to order 10""' the correlation coefficient D{Ee) is commensurable with zero. Hence, to order 10~* the 
imaginary part of the axial coupling constant Aeff is also commensurable with zero, i.e. ImAoff = 0. Nevertheless, 
setting Aoff = RcAoff = —1.2750 and using the relation ImAoff = we may obtain that Im(a^j^) = —0.12 Im(a^^). 
This adds an additional phase shift e'0-S8im(ai^) qj^^ matrix element (Kd)cff- 

The axial coupling constant Xa and the correlation coefficient oq may be also determined by measuring the proton- 
energy spectrum Ea.p3[) (see [40|). The proton-energy spectrum, taking into account the contributions of interactions 
beyond the SM, is 

«(T,) = (l - b,Q)J {.,:(T,) + «o(l + ^ ..(T,)}, (62) 
where the functions gi(Tp) and .92 (7p) SuTC defined by the integrals 



51 (Tp) = "'"^ (1 + 6^ ^) (1 + ^ 9n{Ee)) Ci{Ee,Tp) F{E„Z ^ 1) E, dE„ 

32(Tp)= / "^^'^^ (l + -g.a{E,) + -fn{ES)C2{E,,Tp)F{E,,Z=l)E,dE,, (63) 
where the limits of integration (i?e)max/min ^-re given in Appendix I. 

Summary. We have analysed the sensitivity of the electron-proton energy distribution a{Ee,Tp), the asymmetries 
■^oxp{Ee), Bcxp{Eq), Ccxp and the lifetime of the neutron t„ to contributions of order 10"^ of interactions beyond the 
SM, taken to linear approximation with respect to the Herczeg phenomenological coupling constants of weak lepton- 
nucleon current-current interactions. We have shown that in such an approximation the axial coupling constant Aoff 
and the CKM matrix element (Kid)cff absorb the contributions of the Herczeg left-left and left-right lepton-nucleon 
current-current interactions with the coupling constants a^^ and a'lj^. In this approximation the axial coupling 
constant does not acquire an imaginary part to order 10~^, but the CKM matrix element becomes an additional 
phase e*^™^°J-^+''^«^ Thus, after the measurements of the electron-proton energy spectrum a{Ee,Tp), the proton- 
energy spectrum a{Tp), the asymmetries Acxp{Ee)-, i?oxp(£'e) and Ccxp and the lifetime of the neutron r„ one may 
determine the axial coupling constant Aoff and the real parts of the scalar and tensor coupling constants as and ax, 
defined in Ea. (|i5)) . 

1 + SA^ 

Re(as) +3AoffRe(aT) = —^^bp 

1 + 3A^ 

AeffRe(as)-Re(aT) = t^cst, (64) 
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where in the r.h.s. of the algebraical equations Acff , bp and cst arc the experimental values with their experimental 
errors. 

X. CONCLUSION 

We have analysed the sensitivity of the electron-proton energy distribution a{Ee,Tp), the asymmetries Ac^p{Ee), 
BexpiEe) and Cexp of the correlations between the neutron spin and 3-momenta of the decay electron, antineutrino 
and proton, respectively, for the neutron /3~-decay with a polarised neutron and unpolarised proton and electron to 
contributions of order 10~^ of interactions beyond the SM. For the analysis of contributions of order 10~^ we have 
used the linear approximation for the correlation coefficient with respect to the Herczeg phenomenological coupling 
constants of weak lepton-nucleon current-current interactions. We have shown that in such an approximation the 
Herczeg right-left and right-right lepton-nucleon current-current interactions with the coupling constants a^j^ and 
'^RR Si^^ contributions to the correlation coefficients of the neutron /3~-decay and the lifetime of the neutron. 
Then, the contributions of the Herczeg left-left and left-right lepton-nucleon current-current interactions with the 
coupling constants a^^ and may be absorbed by the axial coupling constant, which we denote as X^g = X — 
Re(a^^ - flii?) + ARe(a^^ + a^^)) and the CKM matrix element (Kd)eff = + Re(a^^ + a^j^)). We have 

shown that the Herczeg coupling constants a^^ and a'lj^ in the effective Hamiltonian of weak interactions may be 
absorbed by the axial coupling constant Aoff = (A — a'lj^ + ci'lj^)/{l + + a^^) and the CKM matrix element 
(Kid)cff = Kid(l + a^L + '^lr)- Using the experimental data on the correlation coefficient D{Ee) we have shown that 
at the level of 10~^ the imaginary part of the axial coupling constant Aoff is equal to zero, whereas the CKM matrix 
element acquires an additional phase e*'™^"^-^-^''^'^-'. This shows that in addition to the background, calculated in 
the SM, the correlation coefficients of the neutron /3~-decay and the lifetime of the neutron, calculated to linear 
approximation with respect to the Herczeg coupling constants of lepton-nucleon current-current interactions, depend 
on the contributions of the scalar and tensor interactions only. This agrees with recent results, obtained in (43l|-j45j. 

We have shown that the contributions of the scalar and tensor interactions beyond the SM are described by the 
Fierz term hp and the coupling constant cst- The Fierz term may be determined from the experimental data on the 
asymmetry Aoxp(i'e) and the electron-proton energy distribution a{Ef^,Tp) (or the proton-energy spectrum a(Tp)). 
The coupling constant cst may be determined from the experimental data on the asymmetry _Boxp(^'e) and the 
proton recoil asymmetry Ccxp- This allows to determine the scalar Re(as) and tensor Re(aT) by solving the system 
of algebraical equations (see Eq.dHH)) 

Re(a5) + 3Acff Re(aT) 
Aoff Rc(a5) - Re(aT) 

The lifetime of the neutron is defined by the background, calculated in the SM, and the Fierz term (see Eg. ([55)) ) 

r„ = (r„)sM(l-6^(^)J, 

where (t„)sm and {rrie/ Ee)sM are calculated in the SM. 

We would like to note that the experimental analysis of the neutron /3~ -decay with a polarised neutron and 
unpolarised proton and electron may be carried out in terms the electron-proton energy distribution a(£'e,Tp), the 
proton-energy spectrum a{Tp), the asymmetries Aoxp(i'e), ^cxp(-£'e), Ccxp and the lifetime of the neutron t„. From 
the fit of the experimental data on these energy distributions, asymmetries and the lifetime we determine three 
parameters, i.e. the axial coupling constant Aetr, the Fierz term bp and the coupling constant cst- In order to 
determine these parameters without correlations between them it suffices to use experimental data, obtained only in 
three from five independent experiments. This means that experimental data, obtained from other two independent 
experiments, should be described well by the parameters (Aefi, &_f, cst), determined from the first three experiments. 
The deviations from the predicted values should be much smaller compared with 10^"', since they may be explained 
only by the contributions of higher powers of the Herczeg coupling constants. 

We have to note that the theoretical analysis of the sensitivity of the electron-proton energy distribution a(i?e, Tp), 
the proton-energy spectrum a{Tp), the asymmetries Aoxp(i?e), BcxpiEe) and Coxp of the neutron /3~-decay and the 
lifetime of the neutron to contributions of order 10~^ of interactions beyond the SM we have carried out above the 
background, calculated within the SM and the account for a complete set of the 1/M corrections, caused by the 
"weak magnetism" and the proton recoil, calculated to next -to-leading order in the large M or the large proton mass 



1 + 3A: 



3Aoff 



CST, 
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expansion, and the radiative corrections of order a/Tr, calculated to leading order in the large M or the large proton 
mass expansion. 

The corrections, caused by the "weak magnetism" and the proton recoil are calculated in analytical agreement with 
the results, obtained by Wilkinson [U and Gudkov et al. [s^]- The radiative corrections to the lifetime of the neutron 
and the correlation coefficients of the neutron -decay with a polarised neutron and unpolarised decay proton and 
electron, calculated in this paper, we have given in terms of two functions {a/n) gn{Ee) and (a/i:) fn{Ee). The 
analytical expressions of these functions are given in Ea. (ID-58p of Appendix D. They are in analytical agreement with 
the radiative corrections, calculated in [l^-js^l and in [31, 32], respectively. We have confirmed Sirlin's assertion that 
an unambiguous definition of the observable radiative corrections to the lifetime of the neutron is fully caused by the 
requirement of gauge invariance of the amplitude of one-virtual photon exchanges of the continuum-state /3~ -decay 
of the neutron [l8|. 

We have improved the theoretical expressions for the asymmetries A^xpiEe), BcxpiEe) and Coxp with respect to the 
expressions, used in [H, 113], [1^, and [42| for the experimental determination of the axial coupling constant A = 
-1.2750(9) and the correlation coefficients A^^"""^^ = -0.11933(34), B^''"''^^ = 0.9802(50) and C(°'^p)o = -0.2377(26), 
respectively. We have added the radiative corrections to the electron asymmetry Ae^piEe) and the 1/M and radiative 
corrections to the antineutrino i?oxp(£'e) and proton Cexp asymmetries. In connection with the experimental analysis 
of contributions of order 10~^ of interactions beyond the SM to the neutron /3~-decay we have calculated the electron- 
proton energy distribution a{Ee,Tp) and the proton-energy spectrum a{Tp) by taking into account the complete set 
of the 1/M corrections, caused by the "weak magnetism" and the proton recoil, and the radiative corrections of order 
a/ir. 

We have also shown that at the present level of the experimental accuracy the lifetime of the neutron is described 
well by the SM with the account for a complete set of the 1/M corrections, caused by the "weak magnetism" 
and the proton recoil, calculated to next-to-leading order in the large proton mass expansion, and the radiative 
corrections of order a/vr, calculated to leading order in the large proton mass expansion. The theoretical value 
of the lifetime of the neutron (r„)sM = 879.6(1.1) s, where the error bars are defined by the error bars of the 
axial coupling constant A = —1.2750(9) and the CKM matrix element Vud — 0.97428(15), agrees well with the 
experimental values Tn^^^ = 878.5(8) s, Tn^^"^ = 88 0.7( 1.8) s and Tn^^"^ = 881.6(2.1) s, measured by Serebrov et al. 
Pilchmaier et al. [76j and Arzumanov et al. |77j . respectively, and the world average values of the neutron 
> ''-^-) = 880.1(1.1) s, T^"-^"-^ = 880.0(9) s and r^"-^"'^ = 881.9(1.3) s, obtained in SSIzi, respectively. 
In order to demonstrate the sensitivity of the lifetime of the neutron, calculated in the SM, to the contributions of 
the radiative and 1/M corrections we propose to rewrite the Fermi integral fn{Eo, Z = 1), given by Ea. ipH)) . in the 
following form 

fn{Eo, Z = l)= j \eo- E,f ^El - ml E, F{E,, Z ^ 1) {(l 



+ ki—gniEe) 

TT 



1 



+k2 TT 



M 1 + 3A2 



(lOA^ - 4(k + 1) A + 2^ Ee - 2A (^A - (k + 1[ 



}dEe, 



(65) 



where the coefficients kj for j = 1 , 2 are equal to kj = or fcj = 1 that means that without kj = and with kj = 1 
corresponding corrections. The numerical values of the lifetime of the neutron for different kj are adduced in Table 
I. It is seen that the most important contributions come from the radiative corrections. 





fci 


k2 


915.3 s 








913.7s 





1 


881.0 s 


1 





879.6 s 


1 


1 



TABLE I: The neutron lifetime, calculated for A = —1.2750, the radiative corrections fci 
k2 = 0, 1, caused by the "weak magnetism" and the proton recoil. 



0, 1 and the 1/M corrections 



For the comparison of the theoretical lifetime of the neutron, defined by Eq.j 
usually used for the measurement of the CKM matrix element Vud [13] (see also [SC 
form M, m 



with the expression, which is 
), we transcribe Ea. ((39| into the 



— = C„]Kd]'(l + 3A2)/(l + RC), 



(66) 



where we have denoted Cn 
radiative corrections [ll, 



= G|m5/27r3 = 1.1614 x 10"'' s'^ and RC = {{a/TT)gn{Ee)) = 0.03886, defining the 
integrated over the phase volume with the account for the proton-electron final-state 
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Coulomb interaction. Then, the phase-space factor /, including the 1/M corrections from the "weak magnetism" 
and the proton recoil, is determined by 

f=^^ f\Eo - E,,f ^El - ml F{E,, z = 1) [l + ^ [(lOA^ - 4(k + 1) A + 2) 

2 

-2A (a - (k + 1)) (^£^0 + ^) }di;e = 1.6894. (67) 

The numerical value agrees well with the value / = 1.6887, calculated in (2^ (see also [lOl)- The factor 1 + RC we 
may represent in the following form 1 + RC = (1 + 5h)(1 + ^r), where Sr = {{a/n) {gn{Ee) — Cwz)) — 0.01505 
is defined by one-photon exchanges only [2^ [2^ and — {a/TT)Cwz ~ 0.02381 is the part of the radiative 
corrections, induced by VF-boson and Z-boson exchanges and QCD corrections [2^, [l^ (see also [13 )• The phase- 
space factor /fl, including the contributions of the radiative corrections, caused by one-photon exchanges only, is 
equal to Jr = /(I + 5r) = 1.71483. It does not contradict the value /a, = 1.71385(34), used in [s^l (see also HH). 

Currently the lifetime of the neutron is proposed to measuring in TU Miinchen within the project PENeLOPE, using 
a superconducting magneto-gravitational trap of ultracold neutrons (UCN) for a precise neutron lifetime measurement 
[stI . In this experiment the UCN are trapped in a multipole field of a fiux density up to 2 T and bound by a 
gravitational force at the top. This makes the extraction and detection of the protons possible and allows a direct 
measurement of neutron decay. A planing accuracy of 0.1 s and better demands high storage times and good knowledge 
of systematic errors, which could result from neutron spin fiip and high energetic UCN which leave the storage volume 
only slowly. Therefore, the neutron spectrum is cleaned by an absorber. The big storage volume of 800 dm"^ and the 
expected high neutron flux of FRMII give more than 10^ neutrons per filling of the storage volume and meet statistical 
demands. Of course, the experimental data on the lifetime of the neutron, which should be obtained within this project 
with a planning accuracy better than 0.1 s, should place new constraints on contributions of interactions beyond the 
SM. 

For the completeness of our analysis we have calculated (see Appendix H) the contributions of the proton recoil 
corrections of order a/A/, caused by the electron-proton Coulomb interaction in the final state of the neutron 
decay. We have shown that these corrections to the lifetime of the neutron and the correlation coefficients are of order 
10~^ — 10~^. This allows to neglect them for the analysis of contributions of order 10^^ of interactions beyond the 
SM. 

We would like to note that we have used the experimental value of the axial coupling constant A = —1.2750(9), 
determined from the experimental data on the electron asymmetry Aexp(£'e) P, H^l- Such an experimental value of 
the axial coupling constant has been obtained with an unprecedented accuracy of about 0.07%. The axial coupling 
constant A = —1.2750(9) agrees well with the axial coupling constants A = —1.2755(13), A = — 1.2759144 5 ^"^^ 
A = -1.2756(30), obtained recently by the PERKEO (PERKEO II) Collaboration (s^ and the UCNA (Ultra- 
cold Neutron Asymmetry) Collaboration [H, [s^ , respectively, the accuracies of which are large compared with the 
accuracy of the axial coupling constant A = —1.2750(9). The lifetimes of the neutron t„ = 879 (2) s, t„ = 879 (6) s and 
T„ = 879 (4) s, calculated for the axial coupling constants A = -1.2755(13), A = -1.2579^4!^:^ and A = -1.2756(30), 

respectively, agree with the experimental data Tn^^'' = 878.5(8) s, Tn^'^'' = 880.7(1.8) s and Tn^'^'' = 881.6(2.1) s, 
measured by Serebrov et al. |l2l |. Pilchmaier et al. [76l | and Arzumanov et al. [77l | , respectively, and the world average 
values of the neutron lifetime t^ ''-^-'^ = 880.1(1.1) s, r^'^ '' ''-^ = 880.0(9) s and r,^'^ ''-''-) = 881.9(1.3) s, obtained in 
Sllllli, respectively. 

The other experimental values of the axial coupling constant A = —1.266(4), A = —1.2594(38) and A = —1.262(5), 
obtained in 0, [sH and [s^, respectively, and cited by Q, lead to the lifetimes of the neutron r„ — 890(5) s, 
T„ = 898(5) s and r„ = 895(7) s. which do not agree with the world average values of the neutron lifetime t^''^'^''^ = 
880.1(1.1) s, tIT''"'^^ = 880.0(9)s and tIT''"'^'^ = 881.9(1.3) s, obtained in i, [zl, [zl, respectively. Moreover, the 
experimental methods, used in [8^ [ssj and j86j for the measurements of the electron asymmetry Aexp(i?e)j has been 
recently criticised in [s^l • As has been pointed out by Mund et al. [s^] , in the experiments [H, |83 and [s^ large 
corrections of about 15% — 30% have to be applied to 1) neutron polarisation, 2) magnetic mirror effects, 3) solid 
angle and 4) background. 

For a long time [l^-flOl (see also [Sll, 113) due to infrared divergences the calculation of the radiative /3~-decay 
of the neutron has been associated with the calculation of the radiative corrections to the neutron /3~ -decay. As 
has been shown already in [l3| . the sum of the rates as well as the electron-energy and angular distributions of the 
continuum-state and radiative /3~~decay modes of the neutron does not suffer from infrared divergences, caused by 
one-virtual photon exchanges in the continuum-state /?~ -decay mode and by the emission of real photons in radiative 
/3~-decay mode of the neutron. 

Nevertheless, the radiative /3~-decay of the neutron n— >p + e^ + Pe + 7 may be treated as a physical process, which 
may be observed separately from the neutron /3~-decay n — > p + e~ 4- z^e- For the first time, the theoretical analysis 
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of the radiative /3 -decay of the neutron n ^ p + e + ly^ + J sls a, physical observable process has been carried 
out in [4^, [5^. First reliable experimental data on the branching ratio of the radiative /3~-decay of the neutron 
gj^(oxp) ^ 3 13(35) X 10^3^ measured for the photon energy region Wmin = 15keV < uj < Wmax = 340 keV, have 
been reported by Nico et al. [i^ . Then this result has been updated by Cooper et al. [i^, [3 1 have obtained 
gj^(exp) ^ 3 09(32) X 10"^. These experimental values agree well with the theoretical value BR^-^ = 2.85 x 10 , 

calculated by Gardner within HB^PT for the same photon energy region [i^, [13, [s^l • In Appendix B we have carried 
out the calculation of the rate, the electron-photon energy and photon-energy spectra and angular distributions of 
the radiative /3~-decay of the neutron with a polarised neutron and unpolarised decay particles. Our results for the 
branching ratios BR^-^ = 2.87 x 10"'^ and BR^-^ = 4.45 x 10"'^, calculated for the photon energy regions Wmin = 
15 keV < cj < 350 keV and Wmin ~ 5 keV < w < Eo — rrie, respectively, agree well with the results BR^-^ = 2.85 x 10~^ 
and BR^-^ ^ 4.41 x lO^'^, obtained by Gardner [i^ and Bernard et al. [53], respectively. Within one standard 
deviation the branching ratio BR^-^ = 2.87 x 10~^ agrees also with the experimental data [46l [47^. 

The rate of the radiative /3~-decay of the neutron in dependence of a photon polarisation has been calculated 
in [50I ]. We argue that the more precise theoretical and experimental analysis of the energy spectra and angular 
distributions of the radiative /3~-decay of the neutron in dependence of the polarisations of the neutron and photon 
should be of great deal of importance for a test of the SM. We are planning to perform such a theoretical analysis in 
our forthcoming publication. 

The theoretical analysis of the sensitivity of 1) the electron-proton energy distribution and the proton-energy 
spectrum, 2) the electron, antineutrino and proton asymmetries of correlations between the neutron spin and the 
3-momenta of the decay particles and 3) the lifetime of the neutron is carried out according to the experimental 
program on the contract I534-N20 PERC, the theoretical program on the contract I689-N16, supported both by the 
Austrian "Fonds zur Forderung der Wissenschaftlichen Forschung" (FWF), and the experimental program of the 
experimental group of the Petersburg Nuclear Physics Institute (PNPI), headed by A. P. Serebrov, on the contract 
No. 11-02-91000 -ANF_a, supported by the Russian Foundation for Basic Research. 
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Appendix A: Amplitude of continuum-state /? —decay of neutron with "weak magnetism" and proton recoil 

corrections to order 1/AI 

The amplitude of the continuum-state /3~-decay of the neutron we rewrite as follows 

G 

M{n ^ pe~ve) = VudMp-, (A-1) 
where Mp- — [up0^iUn]['Se7^(l — 7^)wp] and the matrix takes the form 

= 7^(1 + A7^) + I ^ a^^ikp - knY. (A-2) 
In therms of the time and space components of the matrix = {O^, — O) the amplitude A^^- is defined by 

Mp- = [upO%„] [?2e7°(l - 7^)wp] - [upOun] ■ [?2e7(l - 7^)wp]- (A-3) 
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The time O'^ and spacial O components of the matrix we determine to first in the large M expansion. They read 

1 A + (a ■ L) . 

0^-{ ^ 1 (A-4) 



and 



(A-5) 



where we have kept the terms of order 1 /M only. For the calculation of the amplitude of the /? -decay of the neutron 
we use the Dirac bispinorial wave functions of the neutron and the proton 

u„(d,cr„) = \/2mn( ^" j , Up{kp, <Tp) = y/Ep + nip a ■ kp , (A-6) 

\ Ep + rup^'' J 

where the Pauli spinorial wave functions ipn and (pp depend on the polarisations cr„ and tXp, respectively. The matrix 
elements [upO'^Un] and [upOun] are equal to 

[upO^Un] = ^2mn{Ep + mp) ^^ylifn] + ^ [fli^ ■ kp)ipn]^ (A-7) 

and 

[upOun] = 2mn{Ep + nip) ^^X[^pla (fn] + i ^ [fili^ x kp)ipn] + ^ [fili^ ' fcp)^V«]}, (A-8) 

where in curly brackets we have kept the contributions of the terms of order 1/M only. Using the relation [a ■ kp)a = 
kp + i {a X kp) we rewrite the r.h.s. of Ea. (|A-8p as follows 



[upOu„] ^2mn{Ep + nip) |A[(y9],a^„] + i [ipl{a x /cp)cp„] + ^ bj,"^™]!- (A-9) 

Thus, the amplitude M^- is given by 



M 



^- = ^2nin{Ep + nip)^^[tplip„][ue-/°{l - 7^)wp] - \[ipla ip„] ■ K7 (1 - 7^)up] + ^ [tpl{a ■ kp)ipn 



X [Se7°(l - 7')«.] - * ^ X 4)^„] . [S,7 (1 - j')v,] - ^ [y,t^„] . [fie7 (1 - 7^)^;,]}. (A-10) 

For the transformation of the last term we use the following identity 



based on the Dirac equation for the electron and antineutrino. Substituting Ea. (|A-lip into Eq. (jA-10p we obtain 
^p- = y^2TO„(£'p + mp) + ['Ue7"(l - 7^)wp] - X[tpla (p,,] ■ [uel - 7^)wp] + ^ [tflia ■ kp)ip„] 

X [^Ie7"(l - - i ^ X kp)ip„] ■ [u,j{l - j')v^] - ^[^l^^][u,{l - j')v^]}. (A-12) 



The next step of the calculation is to expand the normalisation factor y/Ep + nip of the bispinorial wave function of 
the proton. This gives 

^2nin{Ep + nip) = 2m„ (l - ^) , (A-13) 
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where we have kept the next-to-leading terms in the large M expansion. 

Thus the amplitude -Mp--, calculated to next -to-leading order order in the large M expansion, is 

^13- = 2m„ |[v3];v3„][ue7"(l - 7^)up] - A[(^t<T(p„] • [?Ie7(l - 7^)wp] + ^ [tflia ■ kp)ip„] 

X [Ue7°(l - 7^)^'P] - * Hi^ X kp)ipn] ■ [Ue7 (1 ' 7^)«P] - ^H^n][Ue{l ~ 7^)«p]}, (A-14) 

where we have denoted A = A(l — Eq/2M). The hermitian conjugate amplitude takes the form 

■^Ij- " 2m„ |[(/3T(y9p][up7"(l - 7^)Me] - ~X*[iplaipp] ■ [wpf (1 - 7'^)lte] + ^ • fcp)(Pp] 

X [up70(l-75)u,] +i!l±-l [<^t(a X • [5^7(1-75^] - Z^[<^t<^J[^;p(l+^5)w^]|. (A-15) 

Substituting Ea. (|X^ into Ea. (fFT|) we arrive at the amplitude of the continuum-state /3 -decay of the neutron, 
taking into account the contributions of the "weak magnetism" and the proton recoil, calculated to next-to-leading 
order in the 1/M expansion. 

The electron-energy and angular distribution is proportional to ^ X^poi l-'^/?" where we sum over all polarisations 
of the interacting particles. Recall that the antineutrino is polarised in the direction parallel to its 3-momentum. The 
quantity i X^poi 1^ equal to 

\M-,-? 

J2 4^2 = t^'il + ^" • o'Ml^e + me)7°fc7°(l - 7')} - Atr{(l + C„ ■a)a}- tr{(fce + me)^k^°{l - 7^)} 

pol " 

-A*tr{(l + f„ • a )a} • tT{{ke + m,h"k^{l - 7')} + |Aptr{(l + 6, • a)aV^}tr{(fc, + m,)j^y{l - j'')} 
tr{(l + 6. • <?)}tr{(fc, + m,h°k{l + 7^)} - ^ tr{(l + Cn • ? )}tr{(fc, + m,)k^"{l 7')} 

+A^tr{(l + en • 0^)0^} • tr{(fce + me)7 fc(l + 7')} + ^^''{(^ + • a )<?} • tr{(fce + me)fc7 (1 - 7')} 
+ ^ tr{(l + U kp)} tr{(fce + me)7°fc7°(l " 7')} 

+ ^ tr{(l + U ■<^){S- kp)} tr{(fce + me)7°fc7°(l - 7')} 
-^tr{(l kp)a} ■ tr{(fc, + m,)^kf{l ~ ^^)} 

-^tr{(l + e; ■cj)cj{<T ■ kp)} ■ tr{(fce + ™e)7°fc7(l - 7')} 

+ • '?)(^ X 4)} ■ tr{(fce + me)7fc7°(l - 7')} 

+ 6. • ^)('? X kp)} ■ tr{(fce + me)7°A:7 (1 - 7^)} 
A*^tr{(l + C ■ d)o={d X kpY} tr{(fce + me)7' W ~ 7')} 
-*A^l^tr{(l + C„ •a)(ax fcp)V}tr{(A;e+m,)7^'fc7^(l-75)}. (A-16) 
Since in our analysis the axial coupling constant A is real, below we set A* = A. Calculating the traces with the real 
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axial coupling constant we obtain the following result 



E 

pol 



(1 + 3A2) EeE + (1 - A") (ke-k)- 2A(1 + A) • ke)E - 2A(1 - A) (C„ • k 



+ — [ - miE + Am^(e„ • fc ) - A (e„ • fee + C„ • fc ) {E,E + k,,-k) + X^ \^{Ee - e„ • ke){h, ■ k + E') 

+ {E + e„ • k){kl + fee • k)] - + 1) [(fn • fce)(fce ■ k + E"") - (fn • fc )(fce' + ' fc ) 

-2(k + 1) A[£'(fc2 + ke-k)- Ee{ke -k + E"^)] + 2(k + 1) A(f„ • fee + fn • k)EeE 



(A-17) 



Taking into account the contribution of the phase volume Eq.(|3]) (see also Eq. (|A-23[) ') and keeping only the terms of 
order 1/AI we have 

\M 

ike,k)J2 3^^2 = (1 + + (1 - A^) (fee ' fc) - 2A(1 + A) (!*„ • fce)S - 2A(1 - A) (!*„ • k)E, 



pol 



"iX^EeE + A2(fce • fc ) + A(2A + 1) • ke)E - A(2A - 1) • /c )Se E^ + 3(1 + 3X^)E,E 



+3(1 - A2) (fee • fc) - 6A(1 + A)(e„ • 4)S - 6A(1 - A) ■ k)E, 

[k-kf ■ ■ 



-3(1 + 3A2)(fce ■k)Ee -3(1 - A^) 



^ +6A(l + A)(ei-fce)(fce-A:') + 6A(l-A)(e„-fc)(fce-A:)^ - m^S + A m2(^„ • A: ) 

-A (e„ • fee + e„ • fc ) (SeS + fee • ) + A2 ^{E, - e„ • fee) ' k + E^) + {E + ' k){kl + k,-k) 

-{k + 1) ■ ke){ke ■ k + E^) - ($„ ' k){kl + 4. • fc)] - 2 (k + l)x\E{kl + k^-k)- E^iK ■ k + E^) 
+2{k + l)A(e, • fee + in ■k)E,E-{K + l)\ ^{i, ■ k,){k, ■k + E^) + ■ k){kl + 4 • ^ )] }■ 
The function (/sg, A; ) is defined by the integral over the antineutrino energy E 



$^-(fce,fc)= / 5{f{E))^ 



n E^dE 



Ep (Eq — Ee 



(A-18) 



(A-19) 



where the function f{E) is f{E) ~ m„ — Ep — Ee — E and Ep = ym^ + [k^ + fc)^ is the proton energy after the 

integration over the 3-momentum of the proton, giving kp ~ —k^ — k. Using the properties of the (5-function the 
result of the integration over E is equal to 



m„ E^ 



Ep {Eq — £'e)2 



df{E) 



dE 



(A-20) 



E = Er 



where Er is the root of the equation f{Er) — 0. To next-to-leading order in the large M expansion the root Er is 
equal to 



Er = (-Bo -E,){1 + ^{E,- k, cos^e^) ) , 



where cost^ep = k^ ■ k/k^E. Since 

TO„ 1 



E„ 



df{E) 



dE 



E = Er 



E^— n 5— = 1 + T7(-^e - keCOS-dev), 

TO„ — iie + Kg cos Vei) M 



(A-21) 



(A-22) 



where we have kept the terms of order 1/M only, using Ea. (jA-21|) and Ea. (|A-22|) for the function $^-(fce,fc) we 
obtain the following expression 



$^-(fce,fc) = 1 + -J^iEe - keCOS-deu) = 1 



M 



E, 



E 



(A-23) 
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We have used this expression for the calculation of Ea. (|A-18p . which defines the correlation coefficients of the 
continuum-state /3~-decay of the neutron, calculated to next-to-leading order in the large M expansion (see Eq.®). 

Appendix B: Radiative /? — decay of neutron 

In this Appendix we calculate the amplitude, the rate, the photon-electron and photon energy and angular distri- 
butions of the radiative /3~-decay n — >■ p -f e~ -t- j/g + 7 of the neutron. The Hamilton operator of weak interactions 
is defined by Eq.([T]), whereas the Hamilton operator of electromagnetic interactions is 

Uemix) = e ['4>p{x)-f'^'4'p{x) - ■0e(x)7''Ve(a;)] Af^{x), (B-1) 

where e is the electric charge of the proton and A^{x) = (0, —A{x)) is the electromagnetic vector potential, taken in 

the Coulomb gauge divA{x) = [ssj . 

For the calculation of the amplitude of the radiative /3~-decay n ^ p + e~ + + j we take into account the 
contributions of the intermediate proton and electron states and drop the "weak magnetism" and proton recoil 
corrections. This gives 

G 1 

M{n-^ pe-Del) = e Kd [up7^(l + A7^)w„] ^^-^ [ue(2e* • fee -f £*g) 7/^(1 - 7^)wp] 

- e^Vud [upi^e* ■ kp + e*q)Y\l + Xl^)un] [uel^il ~ l^H], (B-2) 

where e ~ {0,s) and q = (w,(f) are the polarisation vector and 4-momentum of the photon, obeying the constraint 
e-q = 0. 

We calculate the amplitude of the radiative /3~-decay of the neutron to leading order in the large proton mass (or 
large M) expansion. This agrees well with 1) the iieglect of the contributions of order (a/n) (Eq/M) ^ 10~^ (see 
the discussion below Ea. ipO)) ) and 2) the assertion [13, HI], that the electron-photon energy spectrum, the photon 
polarisation observables and the rate of the radiative /3~~decay of the neutron are dominated by the electron emission 
of photons. In such an approximation the amplitude of the radiative /3~-decay takes the form 

M{n-^pe Vel) ^ e—^Vud ^ , (B-3) 

V2 ^ Ee~n-ke 

where n = q/uj. The amplitude M.^-^ and its hermitian conjugate arc determined by 

Mp-^ = [^l^n][UeQf{l-l^)vA-\ypSVn] ■ [u,Q^{l~l')vA. (B-4) 

and 

M^^^ = [^t v'p][z;,7°Q (1 - 7')".] - Ab,\a^p] • [v,iQ{l ~ 7')^, (B-5) 

where Q = 2{e* ■ k^) -f e*q and Q = -f°Q''j° = 2(e • fc^) + qe. 

The squared absolute value of the amplitude Ea. ()B-4[ ). summed up over the polarisations of the proton and the 
electron in the final state accounting for the polarisation of the neutron is given by 

E = + • '^)} tr{fceQ7°^7"Q(l - 7')} - Atr{(l + • a )(? } • tr{^eQ7°fc7 0(1 " 7')} 

pol. 

-Atr{(l + e„ tr{kQih°Q{l - 7')} + A'tr{(l + • a )a"a"}tr{fce07"fc7" Q(l " 7')}- (B-6) 

Calculating the traces over the nuclcon degrees of freedom and using the properties of the Dirac matrices 

777 =75 -75 +7 5 Ifil , (B-7) 

where £°"^t^^ is the Levi-Civita tensor defined by e^i^s _ ^ ^^^^ eau/ip = —e°"^^^ [s^, we transcribe the r.h.s. of 
Ea. (|B-6p into the form 



Y.\M,-^\' = 2i;[((l + 3A2)-2A(l-A)^)trUeQ7''Q(l-7')} 

pol. 

+ ((1 - A^) -| - 2A(1 + A) C ) • tr{A:,Q7g(l - 7^)}] . (B-8) 
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The traces in Ea. (|B-8[) are equal to 

^ trikQYVil ~ 7')} = (£* • fce)(e • fee) (fee + 9)^ " ^ (£* • e) ((^e • g) - ^ fce^) 

- i (e* • fc,) + (e • fee) (A:e • g + ^ g') + i z e^"^'^ ((e* • fee) £«-(£• fee) £* ) 9^ifce. 

+ ^ « [f <1P-\<1^ ff^') e""^'' 4 £/3 fee., (B-9) 

where = for a real transverse photon. Summing up over the photon polarisations we obtain the following 
photon-electron energy and angular distribution of the radiative /3~ -decay of the neutron 
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For the unpolarised neutron Ea. (jB-10|) coincides with the spectrum, adduced in [l^l (see Eq.Q of Ref.[47j and a 
comment in (89|). 

After the integration over the directions of the photon momentum the photon-electron energy and angular distri- 
bution takes the form 
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Integrating over the phase volume of the final state we obtain the rate of the radiative 13 -decay of the neutron 
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The lowest photon energy Wmin may be treated as the photon energy threshold of the detector. 

For the photon energy interval Wmin < w < w^ax the rate of the radiative /3~-decay of the neutron reads 
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For the region of photon energies 15 keV < uj < 340 keV, used in the experiments by [46|, |47[ , the branching ratio is 
equal to BR^- = 2.87 x lO"^. This result agrees well with the experimental values BRo- = 3.13(35) x 10"^ and 



BR^-^ = 3.09(32) X 10"^ and the result BR^-^ = 2.85 x IQ-^ calculated by Gardner j46 

For the comparison with the analysis of the radiative /3~-decay, carried out by Bernard et al. [sof . we calculate 
the branching ratio for the photon energy region 5keV < uj < Eq — me- The result BR,,- = 4.45 x 10"^ is in a 
agreement with BR^-^ = 4.41 x 10"^^ calculated in [501 ■ 
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For the electron-energy and angular distribution of the radiative /3 -decay of the neutron we obtain the following 
expression 
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where the functions g^^l {Ee,t^min) and g^^j' (i?e,w,nin) arc defined by 
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We would like to note that as we show below the functions g^-^l (i?e,Wniin) and g^-^l (i?e,Wmin), calculated by means 
of the infrared cut-off regularisation Wmin, have energy dependencies different in comparison with the functions, 
calculated in the FPM regularisation (see also f25|'). 

The discrepancy between the electron-energy and angular distributions, obtained for a real photon emission with 
the infrared cut-off and FPM regularisation, respectively, we discuss first in terms of the logarithmically divergent 
integral, which defines the infrared divergent contribution to the amplitude of the radiative /3~-decay of the neutron. 
This integral is 
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where go = v^^+T? is an energy of a photon with mass n, (3 — ke/E^ and v = q/qo are the velocities of the electron 
and massive photon, respectively. The region of the integration in Ea. (|B-16l) is restricted hy qi < q < q2, where 

For the first time the integral Eq. (jB-16|) has been discussed by Kinoshita and Sirlin in It is obvious that 

replacing the lower and upper limits of the integral by Wmin and Eo — E^, respectively, and setting fi zero in the 
integrand the function J {13), given by 
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defines the infrared divergent part of the functions g^^l {Ee,ujinin) and g^^l (i?e,i^min) in Eq. (|B-15|) . 

Now let us calculate the integral Ea. (jB-16p for ^ ^ 0. Following Kinoshita and Sirlin [l^, we may rewrite the 
integral Ea. (|B-16|) as follows 



/32 

J{l3;q2,qi) = y 



+1 



dx 



91 



«^ dqq^ 1 - v'^x^ 
ql (1 - l3vxY ' 



Making a change of variables g — v, proposed by Kinoshita and Sirlin [l^, we arrive at the integral 
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Integrating over x 
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and over v we obtain 
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where is the Spence function, defined by [isl, |9( 
For qi ^ /i the function J(/3) reduces to the form 
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One may see that non-trivial finite and energy-dependent terms appear in the function J {(3), given by Ea. (jB-2ip . for 
gi = and the integration over the region ^ ^ ^ ^2- To show tliis wc set qi — Riid. Q2 — Q^max 

in Ea. (jB-18p and 

arrive at the expression 
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Due to the relation between Spence's functions [19|, |2l|, l9C 



L 



2(3 
1 + /3 



2/3 
1-/3 



2l{ 



2/3 
Vl + /^ 



we derive the following expression for J(/3) 
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we reduce Ea. (jB-26p to the form, obtained by Kinoshita and Sirlin (see Ea. (jC-4p of p^). 



As a result, the functions g^^l {E(.,n) and g'^j^l (i?e,A*), having the form 
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agree fully with the functions, obtained in [23 (see Eqs.(12) and (13) of Ref.[2^). We note that the same result, 
given by Ea. (|B;j6| . for the logarithmically divergent integral Ea. (jB-16p may be obtained within the dimensional 
regularisation [2l|. We will use the functions Ea. (jB-28|) for the calculation of the radiative corrections to the rate and 
the correlation coefficients of the /3~-decay of the neutron (see Appendix D). 
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FIG. 3: Feynman diagrams of radiative corrections to the continuum-state P -decay of the neutron, caused by one-virtual 
photon exchanges. 

Appendix C: Analysis of infrared divergences of radiative corrections to continuum-state /?^— decay of neutron 

The radiative corrections, caused by one-virtual photon exchanges, are described by three irreducible diagrams. 
They are shown in Fig. 3. The diagrams in Figs. 3a and 3b define the self-energy corrections to the masses and wave 
fimctions of the proton and electron respectively. As has been shown in [l^-d^j the self-energy corrections 

to the masses my be removed by the mass renormalisation, whereas the contributions of the self-energy diagrams 
to the wave functions cannot be removed fully by renormalisation of the wave functions of the proton and electron 
and give some observable terms. Wc calculate the contributions of the self-energy diagrams in Appendix D. In this 
Appendix we analyse the contribution of the vertex diagram in Fig. 3c. 

The contribution of the vertex diagram to the amplitude of the continuum-state /?^-decay of the neutron may be 
written as (see Appendix D) 

M^''\n^pe-v,)^e^^V^dM^'^\ (C-1) 



The amplitude M.^'^^ is defined by 
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(C-2) 



where we have kept only the term, which suffers from the infrared divergences. The detailed calculation of the vertex 
diagram in Fig. 3c is given in Appendix D. 

For the calculation of the integral over the 4-momentum q we follow the standard procedure, using Feynman's 
parametrisation of the integrals [3 (see also [9l[). We perform the calculation of the 4-momentum integral in 
Ea. (|C-2p by using the finite-photon mass regularisation p^-fTsj. 



A. Finite— photon mass regularisation 

Below we calculate the integral in Eq. (|C-2p by using the finite-photon mass (FPM) regularisation of the infrared 
divergences. Applying Feynman's parametrisation j9l| we transform the integral in Eq. (jC-2p into the form 

1 /" ^ f ^ f di^ 1 
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4-momcntum integrals [18|,|91[. For the derivation of Ea. (jC-3|) we have used the relation |91l 



where p{x) = kf.x — kp{l — x) and fj, is an infinitesimal photon mass, introduced for the infrared regularisation of the 
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where A^ = - 2ke ■ Q + iO, Ap ^ + 2kp ■ Q + iO, A = - 2p{x) _Q + «0 and B ^ - n'^ + iO, respectively. 
Making a shift of variables q — p{x)y — > g, a Wick rotation qo — iq4 [9l| and integrating over y we arrive at the 
expression [isl IToj 
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where p2(x) — vi^x^ +ina^{l — x)"^ — 2memp"fx{l — x) with 7 — l/\/l — /S^- The r.h.s. of Ea. (|C-5|) mav be represented 
in the form 
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with p = irLp/me- Following then [isl [l9l|. we make a change of variables a — x = b coth(y9. This gives 
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For the calculation of the last integral in Ea. (jC-8|) we make a change of variables i^s = — ^ £nt. This gives 
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where the last two terms are the Spence functions, defined by Ea. (jB-22p . 

Keeping the leading order contributions in the large k expansion we arrive at the expression 
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In terms of the electron velocity /3 it reads 
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where = me/\/l — /3^. For further transformation of the r.h.s. of Ea. (jC-12[) we use the following relation for the 
Spence functions [131 



1 + 13) V1 + /3/ \1-I3J \1-I3J 6 



Substituting Ea. (|€^ into Ea. (|€^ we obtain 



J(/?;/.) = ^^WI^)^nf^)+-Vfi±^)-Lf^)|. (C-14) 
^^'^'^ 2,2tt^ EeTUpPX \ fi J V1-/3/ 4 V1-/3/ V1 + /3/J ^ ' 

Thus, the integral under eonsideration is equal to 

/■ 1 
^ ''^^>- J {2tt)H [q^ + iO] [92 _ 2k, -q + iO] [g2 + 2kp ■ q + iO] 

This result agrees with the expression, obtained in [Hi and [H (see Eq.(B21) of Ref.[lil). 



Appendix D: Total contribution of one— virtual photon exchanges to continuum-state /3 —decay of neutron 



The contributions of one-virtual photon exchanges to the continuum-state /3~-decay of the neutron are shown in 
Fig. 3. The correction to the amplitude of the continuum-state /3~-decay of the neutron, caused by one-virtual photon 
exchanges, we represent in the following form 



where the amplitudes mS^^ and M.Ve define the contributions of the self-energy diagrams of the proton and electron 
in Fig. 3a and Fig. 3b, respectively, and A^pl'' is defined by the vertex diagram in Fig. 3c. Following [l^ we calculate 
the diagrams in Fig. 3 in the Feynman gauge. As a result they are determined by the following analytical expressions 
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where {Snip, Z2 ) and {Srrie, Z2 ) are the renormalisation constants of the masses and wave functions of the proton 
and electron, respectively. 

The dependence on the electron energy is defined by the amplitude M.^pJ only. For the calculation of M^pJ we 
reduce it to the form 
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where we have set = 47ra. The numerator of the integrand of A4pJ we transcribe into the form 
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where we have used the identity For the transformation of the product 
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Substituting Eq. (|D-6p into Ea. (|D-4p we arrive at the following expression of the numerator of the integrand of AA'^ 
in Eq.jEll 
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For the subsequent calculations it is convenient to represent the amplitude Ai^pe as a sum of two contributions 
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and 
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respectively. Now let us proceed to calculate A^J^e; given by Ea. ()D-9[) . Since, the last term in Ea. ()D-9[) is calculated 
in Appendix C, so we should calculate the first three terms only. Using the Pauli-Villars regularisation for the 
ultra-violet divergent integrals (ssj 
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g2 + iO q^ + iO g2 _ ^2 + iO g2 + iO g2 _ A2 + iO ' 
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where A is an ultraviolet cut-off, then Feynman's unification of the denominators, the shift of the virtual 4-momentum 
and the Wick rotation [3^1 for the first three integrals in Ea. (|D-9[) we obtain the following expressions 
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For the calculation of the third integral we have replaced p'^{x) = {keX — kp{l — x))^ by TOp(l — x)'^ 
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For the calculation of the amplitudes J^^ee and A4^) we have to perform the following standard transformations 
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This results in the relation 
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-{q^ + m'^x'^Y {q"^ + m'^x'^Y 
where to = TOe and to = to^ for the electron and proton self-energy contributions, respectively. The term, which is 

.(7) 
Ipp 



not proportional to (to — fc), can be removed by the mass renormalisation. So the amplitudes A^i? and A4^J are 



defined by the contributions of the terms, proportional to (toe — k^) and (rup — kp), respectively. This gives 



= K7^(l + A7^K][Se7,(l-7>.](£/(me) + ^^), 

Hp = K7^(l + A7'K][«e7M(l"7'K](^/K) + ^^), 
where I{m) is given by 
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for TO = TOg and to^, respectively. For the ultra-violet and infrared regularisation we use the following expressions for 
the photon Green function 
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The regularised function /(to) takes the form 
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(^2 + to2.t2)2 (g2 _|_ jji2j.2 _|_ j\2('2 _ 2;))2 

8to2x(1-x2) 8m2x(l-x2) 



(92 + m2x2 + /i2(^l _ 2:))3 {q'2 +m'^x'^ + A^{1 - x)Yi' 18) 

The calculation of the integrals in Ea. (jD-18|) runs as follows. For the first two integrals, which are ultra-violet 
divergent, we obtain the expression 
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The third and fourth integrals in Ea. (jD-18|) are ultra-violet convergent. Moreover in the limit A — > 00 the fourth 
integral vanishes. The third integral is infrared divergent and its calculation gives 
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Thus, the function /(m) takes the form 



A 



I{m) = -2en[ — - Mni — - -. 
Summing up the contributions of the amplitudes Ai^ee , Mpp and A^pl' we obtain the following expression 
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After the cancellation of some terms we arrive at the expression 
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Then, we transcribe the r.h.s. of Ea. ()D-23|) into the following symmetric form 
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Defining the constants of the renormalisation of the wave functions of the electron and proton in the standard form 
as [si 
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we reduce the r.h.s. of Ea. (|D-24|) to the form 

■Mi]^ + M'p'J + M'p'p^ = [Sp7ni + A7')«n] [upi,{l - 7^)1;, 

3 „ /m„\ 1 / Ai \ r 1 „ fl + P 
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In the non-relativistic approximation for the proton Ea. (jD-26|) reads 

'^'^ ^ 27r L V TO,, / 2 V nip / 4 V to^ 
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Now let us take into account the contribution of the first two terms in 5Mpe\ given by Eq. (jD-10p . The first term in 
Ea. (|D-10() is equal to 

<5(i)A^(l) = -(A + 1) ^ ^7^(1 + ^'W][up^^{l - j')v, 
= -(A+l)" 
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TT^i g2 _ 2fcp • q + iO + 2fce • g + iO 
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where we have used Eq. (|D-12p . In the non — relativistic approximation for the proton it reads 

^ -2m„ (A + 1) £ [2^"(;^) + l] {K^n] [uel'il - 7')^'.] - U'^ 'Pn] ■ [uelil ~ l')v,]} . (D-29) 
The calculation of the second term in Eq. (|D-10p runs as follows 
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After the shift of variables, the integration over the 4-dimensional solid angle and the Wick rotation we arrive at the 
expression 
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After the integration over q the term IS 
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(D-32) 

Keeping the leading terms in the large mp expansion, i.e. replacing p'^{x) — {k^x — kp{l — .x))2 by p'^{x) — > '71,2(1 — a;)2, 
and integrating over y and x we obtain 



S^^^MlZ^ = {X + l)^[upr{l + l>n][uel^il~l')v,]{gaf3[in[^^ - .gocffo^ 
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In the non-relativistic approximation for the proton Ea. (jD-33|) takes the form 

S^'^MiV = 2"i"£{(A + 1) en[^) + i] } b>„][Se7°(l - 7')^.] 
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The sum of the contributions of S^^'>M^r>l' and S^^'fM^pJ is 



+ + + J]}['^P^</'"^["«7(l-7'K]• 

Now we have to calculate the contribution of the last term S^'^^AdjZ' to the amplitude 6M]Z\ given by 
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Following the procedure expounded above, we reduce the r.h.s of Eq. (|D-36p to the form 



'■^ dx i 
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where p^ix) = {k^x — kp(l — x))'^ = mlx'^ + mp{l — x)^ — 2mempjx{l — x). 

The calculation of the integrals over x in Ea. (jD-37p has been carried out in detail in Appendix C. Using these 
results we obtain 
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where the ellipses denote the contributions of the terms of higher order in the large mp expansion. 
Keeping the leading terms in the large mp expansion, the contribution of S'^^'^M^J is 
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where we have used the Dirac equation for the free electron Ue{ke • 7) = We(^'e7*' — fne ). Thus <5(3)^^^) 

is given by 
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In the non-relativistic approximation for the proton Ea. (jD-40|) reads 

<5(3)A^(1) = 2m„ £ { [^t^„] [Se7°(l - 7')«.] - Abj;a^„] • [i7e7(l - 7')^.]} 
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The second term can be identified with the contributions of the scalar and tensor lepton-nucleon weak interactions. 
In order to show this we use the Hamilton density operator of weak interactions, taken in the following form [lo| 

nwi:^) - ^Kd{fe(x)7,,(l + A7')V'n(x)] [4(^)7^(1 -7')V'..(^)] 

+ gs [^p(x)^„(x)][V;e(x)(l - 7')V'..(a^)] + ^ 5t fe(x)a^.7'V'n(x)] - 7')^-.. (x)]}, (D-42) 

where gs and gx are the constants of scalar and tensor weak interactions and tr^^ = — 71^7^) is the Dirac 

matrix. 

In the rest frame of the neutron and in the non-relativistic approximation for the proton the amplitude of the 
continuum-state /3~ -decay of the neutron takes the form [loj 



M(n ->p + e + i>e) = "^"^"^ Kd{ [<^^<^„] ["e (1 + gs7°) 7° (1 " 7^)"p] + [vl^ <Pn] ■ [ue(- A + grl") 7(1- 7^)wp } ■ 
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From the comparison Ea. (|D-41[) with Eq. (|D-43p we define the scalar and tensor coupling constants as 
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Using the results, obtained in [l^l, we may define the contribution of the electromagnetic Fierz term [l^] 
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induced by one-virtual photon exchanges. Summing up the contributions of S^^^Mpe with j ~ 1,2,3 we obtain the 
following expression for SM^J 
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Summing up the contributions of Ea. (jD-27p and Ea. (jD-46|) we calculate Aipp +M.ee +M.pe ■ Using this expression we 
get the amplitude of the continuum-state /3~-decay of the neutron, calculated in the non-relativistic approximation 
for the proton and taking into account the contributions of one-virtual photon exchanges. We represent it in the form 

M {n ^ pe-De) = -2m„ ^Kd{ (l + ^ fp- {Ee,fl)) [^l^n] [ue 7°(1 - l^^H] 



-A('-' (1 + ^ fp- {Ee,t,)) [ipla^^] ■ [S,7 (1 - 7')^'.] - £ 9f{E.) H^n] [ue (1 - l')v,] 

+ £ A('^V(i?e)[¥'t?^„] • [^2e7°7(l -7'K]}. 
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Here Gp and g^^ arc the renormalised Fermi and axial coupling constants 
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where dy and cIa are ultra-violet divergent constants. In our calculation of the radiative corrections they are equal 
to 



(rrir,) 4 ^ ^ ^ -2 (mr, 
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The function fp-{Ee, fJ-), given by 

/..(--)-H3-^--(£-)[2^-(l^) 



/?^ll + /?J 4/3^" Vl-/?J + 2/3^"Vl-/?J' 
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contains two terms {3/2)£n(mp/ rrie) and C5 independent of the electron energy E,,. As we show in Appendix E 
the term {3/2)£n{mp/me) is fixed by the KLN theorem. In our calculation the constant C5 reflects an ambiguous 
decomposition of the contribution of one-virtual photon exchanges to the amplitude of the neutron /3~-decay into 
the renormalisation constant dy of the Fermi coupling constant and the radiative corrections to the lifetime of the 
neutron. Nevertheless, following Sirlin [l^ and Abers et al. [l^ one can show that due to a requirement of gauge 
invariance of the observable radiative corrections the value of the constant cs is fixed and is equal to cs = —11/8. 
We show this in Appendix F. 
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The contribution of the VF-boson and Z-boson exchanges and the QCD corrections [2^ we describe by the constant 
Cwz [32] • The numerical value Cwz = 10-249 wc discuss below Ea. (|D-58p . As a result the function fp-{Ee,iJ,) is 
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Making a replacement G^p — >■ Gf and A'-''^ — )- A and taking into account the contributions of the "weak magnetism" 
and the proton recoil the amplitude of the continuum-state /3~ -decay of the neutron takes the form 

M(n ^pe- De) = -2m„ ^Kd{ (l + ^ [Ee, m)) [vlv>n][ue 7°(1 - 7^)«p] 
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where A = A(l — Eq/2A-I) and kp — — fcg — /c is the proton 3-momentum in the rest frame of the neutron. The 
amplitude Eq. (|D-52[) has been used for the calculation of the electron-energy and angular distribution Eq. ^ . 

The radiative corrections to the rate of the continuum-state /3~-decay of the neutron, which we denote as g^- {Ee, 
acquire an additional contribution of the electromagnetic Fierz term 
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where the term, proportional to /3/2, is defined by 
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For the calculation of the contributions of the effective scalar and tensor interactions, induced by one-virtual photon 
exchanges, to the correlation coefficients of the electron-energy and angular distribution of the neutron /3~-decay we 
use the results obtained in [H (see Eq.(28) of Ref.fl^). The corrections to the correlation coefficients from the scalar 
and tensor lepton-nucleon weak interactions with the left-handed neutrinos take the form 

h "^'^ o 9s - 3A.gT me 



5a{Ee) = 0, 
5A{Ee) = 0, 

5B{Ee) = 2 



[gr - Xgs) - '2>^gT rUe 



1 + 3A2 



Ee 



(D-55) 



where we have kept only the linear terms in the scalar and tensor coupling constant expansions. Replacing the scalar 
and tensor coupling constants by their expressions, given in Ea. (jD-44[) . we obtain 
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This gives the foUowing radiative corrections to the correlation coefficients 
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where the terms —gprrie,/ are absorbed by the function gn{Eg) (see Ea. (|D-53p ). The terms of order of 0{1/M) are 
adduced in Eqs.® - (fT5|) . 

The functions gn{Ee) and /ri(-Ee) are defined by 
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where the functions g^^l_^{Ee, n) and g^^l^{Ee, fJ-) are given in Eq. (jB-28|) . 

The gn{Ee) and /„(i?e), multiphed by a/7r, are in analytical agreement with results, obtained in (T^-js^l and 
[31I [3^ . respectively. The constant Cwz-, defined by the contributions of the VF-boson and Z-boson exchanges 
and the QCD corrections [2^, is equal to Cwz = 10.249. This numerical value is obtained from the fit of the 
radiative corrections to the lifetime of the neutron {a/iT){gn{Ee)) = 0.03886(39) [l| and (a/7r)(g„(£'e)) ~ 0.0390(8) 
(28l | , averaged over the phase volume of the neutron decay. 

The radiative corrections {ol/tz) g„ {Ee) and {ol/tt) fn{Ee), weighted with the electron energy spectrum density 



Pp- (E,) = (Eo - E,f ^El - ml E, ({E,) 
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where the functions C{Ee) and F{Ee,Z ~ 1) are given in Eq.© and Eq.(IS]), respectively, and fn{Eo,Z = 1) is the 
Fermi integral Ea. (|40)) . are plotted in Fig. 1. 

Finally we would like to note that having attributed the terms, proportional to (A + 1) to the renormalisation 
constants of the Fermi and axial coupling constants only, we arrive at the radiative corrections, described by the 
function 



g{E,) = g„{E,)+3en(—) + 
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— Cwz — 



= Mn{ — 
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The function g{Ee), multiplied by a/ir, agrees analytically with the result, calculated by Kinoshita and Sirlin [l5| 



Appendix E: Kinoshita-Lee-Nauenberg theorem for radiative corrections of neutron /3 —decay 

As has been shown by Kinoshita and Sirlin [isj . the radiative corrections to the rates of the muon decay fi^ — > 
+ + i?e and of the neutron /3~-decay, taken in the limit rrie — )■ and integrated over the phase volume, does not 
depend on the electron mass. This is so-called the Kinoshita-Lee-Nauenberg (KLN) theorem [3g|. Thus, the radiative 
corrections to the rate of the neutron decay, described by the function gn{Ee), should obey the KLN theorem. This 
means that the result of the integration of the function g„(i5e), taken in the limit nie — >■ 0, over the phase volume 
should not depend on mg. 
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Sirlin's function g{Ee) = p„(i?e) — Cwz, defining the radiative corrections, caused by one-virtual photon exchanges, 
to the hfetime of the neutron, takes the form [13 



g{E,) = g„iE,) -Cwz = \ inC^) - ^ + 2 

2 V TTip / 8 



^lnC-±^\-x\ \ln{ 
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nie 



3 I Eq — Ee 

2 3 



E,_ 



(E-1) 



Taking the hmit Ee ^ me, corresponding to /? — > 1, and introducing the variable x = E^/Eq wc transcribe the 
function g{Ee) into the form 
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Thus, the function g{x,me) is equal to 

3 



g{x,m, 

1 (l-x)2 
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For the derivation Ea. (|E-3p we have used L{1) = —ir'^/G and the approximation 
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According to the KLN theorem [sl] (see also [l^), the function g{x,me), integrated over the phase volume, taken in 
the limit TOp — 0, should not depend on mp. We would like to emphasise that the quadratic terms in^(2E()/me) are 
cancelled in the function g{x, TOp) without integration over the phase volume. 

Now let check the contribution of the linear terms in{2Eo/nie). For this aim we have to integrate the function 
g{x,me) over the phase volume. The integration of the function g{x, TOp) over the phase volume with a dimensionless 
element (1 — x)^x'^dx, obtained at TOp — >■ 0, gives the result independent of Top, since 



2en 



fl-x\ 3 2 1 
(-^)-2 + 3- 



1 (1-^-)^ 
12 a;2 



(l-x)^^^ dx^O. 



(E-5) 



Thus, the term £n(2Eo/me) vanishes. This reproduces the results, obtained in [T3, and confirms the KLN theorem 

We would like to note that the term —3/2 in the integrand of Ea. (|E-5l) . playing an important role for the vanishing 
of the integral, is given by —3/2 = —3 + 3/2. where —3 and +3/2 come from the energy depending part of the function 
g{Ee) and the term {3/2)in(mp/me), respectively. 

Appendix F: Comparison with Sirlin's calculation of radiative corrections |l8|| 



In this Appendix we compare our calculation of the radiative corrections to the continuum-state f3 -decay of the 
neutron, caused by one-virtual photon exchanges, with the calculation, carried out by Sirlin in his well-known paper 

[3- . . . ^ 

According to Sirlin [18| , the radiative corrections to the amplitude of the continuum-state /3 -decay of the neutron 
are defined by Eqs.(9a), (11), (14) and (19) of Ref.[3. Now let us compare our expressions with Sirlin's ones. 

Sirlin's Eq.(9a) corresponds to our amplitude Aipe- In order to show this we rewrite the amplitude A4pe, given by 
Ea. (|D^ . as follows 
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Top 



iO 



1 



kp + q — iO 



(F-1) 
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where Dap{q) is the photon Green function 



in the arbitrary gauge with a gauge parameter ^, and are the products of the Dirac matrices defined by 



(F-2) 



(F-3) 



After some algebraical transformations Ea. (jF-l|) can be reduced to Sirlin's form 



'■pe 



Dap{q) 



where T*^^ is given by 
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q'^ -2kp-q + iO 



(F-5) 



and i?^^ is equal to R^^ = ia^''^q\W'^. As has been pointed out by Sirlin [T^, the tensor T^'^ is obviously transverse. 
g^T''^ = 0. Then, we propose to rewrite the amplitude Eq. (|F-4[) as follows 



(F-6) 



where the indices (SC) and (SLI) mean "Sirlin's Corrections" and "Strong low-energy interactions" . The amplitudes 
A^pe^^ and SM^fe^^^ are equal to 



5_A/((SLI) 



47r J TT^i 
a f d^q 
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The contribution of Sirlin's Eq.(ll) coincides with the part of the amplitude Ai'^eJ ^ defining renormalisation of the 
wave function of the electron by electromagnetic interactions. The amplitude fA^ee , given by Ea. (|D-2[) . we rewrite as 
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where the amplitudes M.^ee^'^ and SJ^fe'"'^' are equal to 
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(SC) 

respectively. For the derivation of Sirlin's term Mee one has to use the definition of the renormalisation constant 
z'^'^ — 1 of the electron wave function [l^ [3^ 
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(F-10) 
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In comparison with Sirlin's expression, we have taken away the operator (fee + '7i,e)/2me, which is equal to unity at 

- — - TJl/Q . 

As has been pointed out by Sirhn, Eq.(14) of Ref. fisj is related to the emission and absorption of a photon by the 
proton. This means that Eq.(14) of Ref.fisj should be a part of the amplitude Ad^^p ■ This allows us to represent the 
amplitude Ai^p in the form 

M%^=Mlf^+5Mf^'^ (F-11) 

with Alpp*"^ and SM^^\ given by 



(p) 
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[upW^UrMu^O.v,] I ^^oMg) ^^!^ oT^l^^jf ' (F-12) 



respectively. As a result, according to Sirlin [l^, the observable radiative corrections to the amplitude of the 
continuum-state /3~-decay of the neutron, caused by one-virtual photon exchange, are defined by the amplitude 

■MI^c ^ = -^r^ + Mil^^ + Mg^\ (F-13) 

which is gauge invariant, i.e. invariant under the gauge transformation of the photon Green function Dap(q) — > 
Dap{q)+c{q'^) QaQfi with an arbitrary function c{q^), and suffers from the infrared divergences only [l8| . The subscript 
RC means "Radiative Corrections" . As we show below the additional contribution, described by the amplitude 

SM^^'^ = SMlf^'^ + SMif'^ + SM^f'^ , (F-14) 

does not depend on the electron energy and should be absorbed by renormalisation of the Fermi coupling constant 
Gf and axial coupling constant A. 

Since the amplitude -M^q is invariant under gauge transformations of the photon Green function Dap{q) — > 
Daf}{q) + c(g2) qaqp, we may calculate it by using the Feynman gauge for the photon Green function. The calculation 

of the amplitude SM^^^^K^ describing the contributions to renormalisation constants of the Fermi and axial coupling 
constants, can be also carried out in the Feynman gauge. 
The amplitude A^pe can be rewritten as follows 
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■kH q^ - 2kp ■ q + iO q^ + 2ke ■ q + iQ ' J -kH q^ + iQq^-2kp-q + iO q^ + 2ke-q + iO 

/ ^ ,2',0g2,2fcp.g + ^0g2 + 2fc!.g + ^0 ["^-"^'^^^p'^,.,]}. (F-15) 



All momentum integrals in Eq. (|F-15|) are calculated in Appendix D by using the Pauli-Villars regularisation of the 
ultra-violet divergent integrals and the FPM regularisation for the infrared divergent contributions. Using the results, 
obtained in Appendix C and Appendix D (see Eas. (jC-15p . (|D-12p and (|D-40|) ). the amplitude Alpl^'' takes the form 
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(SO 

The amplitude Aiee can be transcribed into the form 
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Having integrated over x we arrive at the following expression for the amplitude A4 
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In the Feynman gauge the amphtude A^pp takes the form 
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Using the results, obtained in Appendix D, we define A^pp^-* as follows 
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(F-20) 



After the summation of the contributions the amplitude takes the form 
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Thus, the value of the constant cs is cs = —11/8. It is fixed by the requirement of gauge invariance of the observable 
part of the radiative corrections to the amplitude of the continuum-state ^"-decay of the neutron. 
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Now let us proceed to calculating the contribution of SM.^q^\ The term Ai^J"^'' is defined by 
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where i?^ = 7^(1 + 7^). Using the results, obtained in Appendix D (see Ea. (|D-12p . Ea. (|D-28p . Ea. (|D-33|) and 
Ea. (|ro8)) ) we obtain 
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For A^if'"^'' and A^pp'^^' we calculate the following expressions 
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respectively. Using the definition of the renormalisation constants z'^^ and ^2^'' Eq. ()D-25P we obtain 

TWi^") = 0, 

The sum of the contributions to A^^fc^^ gives one 



(F-26) 



27r 



+ 2;^(^ + l) 



3 to 







f-2tofA) 


- 1 


V nrip / 





1 



(F-27) 



This shows that gauge non-invariant contributions do not depend on the electron energy Ee and the infrared cut-off fi 
and may be fully absorbed by the renormalisation constants of the Fermi and axial coupling constants (see Eq. (|D-48[) 
and Eq.jEMl)). 
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In order to prove the correctness of the term —3/8 in Eq. (jD-59|) . giving the contribution to the observable radiative 

fsc) 

corrections to the Ufetime of the neutron, we propose to sum up the contributions of A^j^^^ and the terms, proportional 
to [upW'^u„][ueOfj,Vi^] from which we denote as A^jfc^^- This gives 
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where the function F(i?e,/i) takes the form 
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The contribution of the radiative corrections to the lifetime of the neutron is defined by the function 



(F-29) 



\im[FiE,,^,)+g^'liE,,^l)] 
A 



2/3 VI -/3 



2(£;o - ^e) 



1 Eq — Ee 

3 



13 J ' 2/3'^"'Vl -/3y L^^ ' " ' ' 12 £:2 



(F-30) 



where the function gj^l{Ee, iJ,)is given by Ea. (|B-28p and describes the contribution of the radiative /3 -decay of the 
neutron. 

The function Ea. (|F-30p reproduces our result, obtained in Appendix D (see Ea. (|D-60|) ). Then, being multiplied by 
(a/Tr), it reproduces also the radiative corrections to the lifetime of the neutron, calculated by Kinoshita and Sirlin 
[isj . This corroborates the correctness of our calculation of A^^c'^ and -M^c^^ ■ 

Thus, the analysis of the radiative corrections to the neutron /3~-decay, carried out in this Appendix, confirms 
Sirlin's assertion that the requirement of gauge invariance of the amplitude of the radiative corrections defines un- 
ambiguously the observable radiative corrections to the lifetime of the neutron. They are independent of the axial 
coupling constant A, i.e. of strong low-energy interactions. The part of radiative corrections depending on the axial 
coupling constant A, i.e. on strong low-energy interactions, is unobservable and absorbed by renormalisation constants 
of the Fermi and axial coupling constants. 



Appendix G: Contribution of weak lepton— nucleon couplings beyond SM to correlation coefficients of 

neutron /3^— decay 



In this Appendix we take into account the contributions of the weak lepton-nucleon interactions beyond the SM 
with left-handed and right-handed neutrinos. For this aim we use the following Hamiltonian of phenomenological 
weak lepton-nucleon interactions 

G f - - - - 

Uwix) = Kd|[V'p(2;)7,xV'n(a;)][Vie(a;)7^(Cy + Cy7^)V'^^(a;)] + [V'p(a;)7^7^Vn(a;)][V'e(a;)7^(C'A + CAf)^yAx)\ 

+ [M^)^n{x)][4^,{x){Cs + Csi'')^.Ax)] + \[M^)^^^ (G-1) 

This is the most general form of the effective low-energy weak interactions, where the coupling constants Ci and Ci 
ioT i = V,A,S and T can be induced by the left-handed and right-handed hadronic and leptonic currents and 
supersymmetric interactions Q as well. They are related to the coupling constants, analogous to those which were 
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introduced by Hcrczeg [6|, as follows 
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Ca 
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_ Ah _ Ah i_/l'' +4'* 
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(G-2) 



where the index h means that the coupling constants are introduced at the hadronic level but not at the quark level 
as it has been done by Herczeg 0]. In addition in comparison with Herczeg Q we have taken away the common 
factor CpVud/V^ and defined the coupling constants a^^ and a^^ as deviations from the coupling constants of the 
SM [o^l- Analogous to Hcrczeg 0, the Hamiltonian of phenomcnological lepton-nucleon weak interactions beyond 
the SM may be written in the following form [92 1 
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■1 - A 
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In order to express the coupling constants Ct and Ct in terms of the coupling constants and we have used 
the relation 0^,^^^ = |e^i/a/30'"'^ [H]. 

The SM is defined by the coupling constants Cs ^ Cs = Ct = Ct = 0, Cy = — Cv = 1 and Ca = — Ca = — A 
and the couphng constants a^, Afl^ and a^^ for = L or R are induced by interactions beyond the SM. 

The amplitude of the continuum-state /3~ -decay of the neutron, calculated with the Hamiltonian of weak interac- 
tions Eq. (|G-ip to leading order in the large proton mass expansion, takes the form 
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The hermitian conjugate amplitude is 
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+ VpIVp\[vu{C*s - C*s'^^)ue\ - yiSipp] ■ [wp7"7(C^ - C|.7^)ue]|. 



(G-5) 



The correlation coefficients of the neutron j3 -decay, expressed in terms of the phenomcnological coupling constants 
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Cj and Cj for j ~ V, A, S and T and calculated to leading order in the large M expansion, are equal to 

^ = |CvP + |Cv|2 + 3|CaP + 3|CaP + |CsP + |CsP + 3|CtP + 3|CtP, 
= \Cv\' + \Cv\' - \Ca\' - \Ca\' + |CtP + |CtP - |CsP - \Cs\', 

= 2Re((CyQ + CvC*s) ~ 3{CaC^ + CaC^)) , 
= 2 Re(2 CaQ - 2CtC:^ - (CyC^ + GvC*a) - (CsCJ + CsCJ)) , 
= -2Re(2CACl + 2CTC; + (Cy(7:^ + C'yC':^) + (C'sC'|, + CsC'|,) 



+2Rc(^(CvCJ + CvC*t) - {CaC*s + CaC*s) + 2 (C^C^, + CaC*t,)) 
CD = 2Im((CyCl + CyCl) + (CsC^ + CsCT))- (G-6) 



The correlation coefficients Ea. (jG-6|) reproduce well the structure of the correlation coefficients, calculated in 

At the neglect of the 1/M and radiative corrections the rate of the neutron /3~-decays, expressed in terms of the 

coupling constants Cj and Cj for j = Y ^ S and T is equal to 



where the Fermi integral fn{EQ,Z = 1) is given by Eg. ipH)) at the neglect of the l/M and radiative corrections. In 
the SM at the neglect of the 1/M and radiative corrections ^ = 2 (1 + 3 A^). 

The correlation coefficients of the neutron /?~ -decay with the contributions of weak interactions beyond the SM 
Ea. (|G-ip . defined as the deviations from the correlation coefficients calculated in the SM with the account for the 
contributions of the "weak magnetism" and the radiative corrections, are equal to 

C{Ee)^CsM{Ee)+(bc+bF'^), 

a{E,) = asuiEe) + \ ((|Cyp - 1) + {\Cv\' - 1) - {\Ca? X') - (|CaP - A^)) - ao [be + bF'^), 

A{E,) = AsuiEe) + Re(2 {CaC*a + A^) - (CyQ + CyQ - 2A)) - Ao (be + bp 

B{E,) = BsuiEe) - Re(2 (C^Q + A^) + (CyQ + CyQ - 2A)) + 

/ — — — — — — \ 771 / 771 

X Rc({CvC^ + CvC^) - {CaC*s + CaC*s) + 2 (CaC^ + CaC^)) - Bq [be + bF^ 

D{E,) = Dsm{E,) + Y^L_Im((CyCi - A) + (CyC^ - A)), (G-8) 

where the correlation coefficients CsM{Ee), asuiEe), AsuiEe) and BsuiEe) arc given by Eas. ([TU| - (|T^ . respectively. 
They are induced by the SM with the contributions of the "weak magnetism" , the proton recoil and the radiative 
corrections. The terms in the correlation coefficients Eq. (jG-8p . depending on the phenomenological coupling constants 
Cj and Cj for j = V, A, S and T, define small corrections to the contributions, calculated in the SM. In Ea. (|G-8[) 
we have dropped the contributions of the squared values and crossing products of the scalar and tensor coupling 
constants. 

The correlation coefficient _DsM(-E'e) has been calculated in (H^-fg^]. For the electron kinetic energies 250 keV < 
Te < 455 keV and the axial coupling constant A = —1,2750 one may estimate that DsuiEe) ^ 10~^ [5^. Then, the 
coefficients be and bp are defined by 

be = \Y^{i\Cv\'-l) + i\Cv\'~l) + 3{\CA\'-X')+3i\CA\'-\')), 

bp = ^-^Re(^{CvC*s + CvC*s)-HCAC^ + CAC^)), (G-9) 

where bp is the Fierz term P, d, 0, Q- -l^o^ ^^.e analysis of the sensitivity of the observables of the neutron /3~ -decay 
to contributions of interactions beyond the SM we define the corrections to the correlation coefficients Ea. ljG-Sp . 
caused by new phenomenological weak interactions Ea. (|G-l|) . in the linear approximation with respect to the Herczeg 
coupling constants a^^, A^^- and af^ for = L,R Eq. (jG-2p (see section HX)) . 
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The rate of the neutron p -decay, corrected by the contributions of interactions beyond the SM, takes the form 

A„ = (A„)sMfl + &c + 6F(^) ), (G-10) 

where (A„)sm is defined by Ea. (|5ni) and {nie/ Ee)sM is the average values, calculated with the electron-energy spectrum 
Ea. (|D-59p . The lifetime of the neutron is equal to 

Tn = {Tn)sM(l-bc-bF('^) ). (G-ll) 
V ^ Ep l SM/ 

The contributions of the terms he and hp are discussed in section HXl We show that the contribution of he to the 
lifetime of the neutron may be removed by the redefinition of the axial coupling constant and the CKM matrix 
element. 



Appendix H: Contribution of proton recoil, caused by electron— proton final state Coulomb interaction 

In this Appendix we calculate the contribution of the proton recoil, caused by the electron-proton Coulomb inter- 
action in the final state of the decay. A velocity of a relative motion of the electron-proton pair is equal to 

v=^-^. (H-1) 
rUp Ee 

The Coulomb corrections, caused by a proton recoil and calculated to order a/M, change the Fermi function Eq.([5]) 
as follows 

, / 7ra Ep TTa En — E„ kp ■ k\ , , 

FiEp^Z^ 1) ^ FiEp^Z = 1) (l - - ^ - ^ ^ (H-2) 

As a result the function C{Ep) acquires the following correction 

CiEp) ^ CiEeMEe) = C{Ep){l - y §(l + 5 «o ^^))- (H-3) 

The function rj{Ep), averaged over the electron energy spectrum pp-{Ep) Ea. (jD-59p and equal to {rj{Ep)) = 1 — 
2.7 X 10~^, defines a correction to the lifetime of the neutron of order 10^^. The correlation coefficients acquire the 
following corrections 

^ ' "3/32 MEp M ' ^ ' " /3 M ' ^ ' " /3 M' 



6K,,{Ep) ^ -A,^-^^^ , 5Q^{Ep)^-B,^-^^^. (H-4) 



In the electron energy region 250 keV < Tp < 455 keV and at A = —1.2750 the obtained corrections are of order 

5a{Ep) 10-^ 5A{Ep) ~ 10-^ 6B{Ep) 10~^ SK„{Ep) - lO"*^ and SQn{Ep) 10"^ Thus, the corrections, 

induced by the proton recoil in the Coulomb electron-proton interaction in the final state of the neutron /3~ -decay 
and calculated to order a/M, are smaller compared with corrections of order 10~^ of interactions beyond the SM in 
the experimentally used electron energy region 250 keV <Tp< 455 keV. 



Appendix I: Electron— proton energy distribution and proton recoil asymmetry of neutron /3 —decay 

The correlation coefficient a{Ep), which analytical expression is given in Ea. ljlip . can be hardly used directly for 
the experimental determination of the correlation coefficient ao due to impossibility to determine experimentally a 
correlation between the electron and antineutrino 3-momenta. For the experimental determination of the correlation 
coefficient ag one needs to measure the correlations of charged particles, i. e. the correlations between 3-momenta 
of the proton and electron. For the experimental analysis of the proton recoil asymmetry we need also to have the 
electron-proton energy and angular distribution, including the correlations between the neutron spin and the proton 
3-momentum. 
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After the integration over the antineutrino S-momentuni the electron-proton energy-momentum and angular dis- 
tribution of the neutron /3~ -decay takes the form 



(I-l) 



where E — \kp + ke\, dQp ~ sin 6pddpd(j)p is the infinitesimal element of the solid angle of the 3-momentum of the 
proton relative to the polarisation vector ^„ of the neutron, i.e. ^„ • fcp = PkpCO&Op with P = < 1. Then, 
d^ep = sinOepdOepdcjjep is the infinitesimal element of the solid angle of the correlations of the electron-proton 3- 
momenta kg ■ kp = kgkpCosOep and ■ ke = Pke{cos Op cos 6(,p + sin^p sin^ej, cos((/)j, — 4>ep))- In the non-relativistic 
approximation for the proton and to next -to-leading order in the large M expansion we replace nin/Ep by 

^^1 + ^. (1-2) 
Ep M ^ ' 

Then, |A^^-p, multiplied by the factor Eq.(|01)), is (see Ea. (jA-17l) ') 

+ i(i?e) + B(i?e) %^ + ^ YtW{ ~ + ^ ■ + ^^^^ 

X • • - (a^ + 1)A + (^ + 1)) • ^j^^^ • }, 

(1-3) 

where A; = —kp — k^, E = \kp-\- kp\ and we have denoted 

1 - A2 M A TT 



a(i?e)=ao(l + ^ r^^) (l + 3/„(ii;e) 



i(i^e)^Ao(l + (l + ^/„(i..)) 



2(1 + A) il/. 

^(^^)-^°(i + 2(T^i^(^-f))- (^-4) 

Substituting Ea. (|I-3p into Ea. (jl-l|) we obtain the electron-proton energy-momentum and angular distribution of the 
neutron /3~-dccay with polarised neutron and unpolarised electron and proton 



dEpdkpd^pd'^^p - (1 + j 1 + «(£^e) ^ + YT3A^ 1^^" ^ ^ " +2(,. + l)Aj-^ 

- (a^ - 2(. + 1)a) + i(i..) + i.(i..) ^ ( - (2. + 1) A (e; ^ 

+ (a^ + + 1)A + (K + 1)) • • - (a^ - (- + 1)A + (. + 1)) (^lilM^^j I 



X (^1 + - .g„(Se)J F(Se, ^ = 1) '5(m„ - Ep - E, - E) k^E^k;, (1-5) 

Now let us analyse the case of unpolarised neutrons. Integrating over the solid angles dflp and dflep we obtain the 
electron-proton energy spectrum 



d^X^-iE^,Tp) _ ,,,,G|,|K.'^ 



a 



^j^^^rj.^ ^M{l + iy) a{E,,Tp)[l + -gr,{E,))F{E,,Z = l)E,, (1-6) 

where Tp = kp/2M is the kinetic energy of the proton. The electron-proton energy distribution a{Ee,Tp) takes the 
form 

a{E,,Tp) = Cl(Se,rp) +a(£;e)C2(Se,rp). (1-7) 
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The functions (^i(Ee,Tp) and C2(-£'e,'7p) are determined by the integrals over cos^ep = ■ kp/k^kp 

^ Kk, I*^{,,L_J_[e„-^, (a= + 2(. + 1)A) - (A= - 2,K + 1)a) ^ } 

X s(jnn — Ep - Ee — \kp + ke\j dcosOep (1-8) 

and 

C2Ee,Tp) = -fcefcp / +ke-kp ^/ - Ep - E, - \kp + fce|) rfcOS^ep- (1-9) 

J-1 Ee\kp + ke\ ^ ' 

The integration over cos Q^p we perform by making a change of variables E ^ | fcp + /sg | . This gives 
Ci(i?e,Tp) = {iJ+l_L-^[(iJo-|i)iJ + i(A2 + 2(K+l)A)(ii;2 + fc^-fcf) 



1Ka^_2(.+ 1)a) (i?^-fc^-fc,2)]}, 
1 1 



UE,,Tp) = -- — {E'^k'p + ki) (I-IO) 
at i? = m„ — Ep — Ee- The energy region of the definition of the distribution a{Ee, Tp) is given by [osj 

0<Tp< (T,)_ = '-^^^ ^ = = 0.751 keV, 

(^'e)min < < (-Eg) max, (I-H) 

where (i^e)min/max are equal to 



_ {m,,-EpTkpY + ml _ {E,T^2MTpf+ml 1( ^1 ^ _ 



2{mn-E„Tkp) 2 (£^0 T 72177;) 2V (E^T y/2M%)^. 

Here we have kept the next-to-leading terms in the large M expansion and used the relation 

m„ - mp = Eo + ^^^^ ^Eo + {Tp)^,^. (1-13) 
To next -to-leading order in the large M expansion the variable E is defined by 

E = m,,~Ep^E,^Ea-E, + {{TpU,^ - Tp). (1-14) 
The functions (i{Ee,Tp) and (2{Ee,Tp), determined to order l/M, are equal to 

(l{Ee,Tp) = (£^0 - Ee) + ((Tp)max ~ + ^ + SA^ ^ ~ 

((i^o - E,f + El -ml- 2MTp) + ^ (a^ + 2(a^ + 1)a) ((£o - E,f - El + ml - 2MTp) 
-\ (}? 2{n + 1)a) ((iJo - E,f - El + - 2Mr,)" , 

CsSe, ^p) = ((^0 - Se)' + - - 2MTp + 2(^0 - i^e)((Tp)„,ax - Tp)) . (1-15) 

Integrating the electron-proton energy spectrum Ea. p-6|) over the electron energy we obtain the proton-energy spec- 
trum 

^^^.M(l + 3A^)%i^a(T,), (1-16) 
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where aiTp) is defined by 



1 E \ 

a{Tp) ^ giiTp) + ao{l + Y^7j2 Jij) 92iTp) 



(1-17) 



with the functions gi{Tp) and g2(Tp), given by the integrals 

r-(i5e)„ax 



9iiTp) 



92{Tp) = 



Ci(-Be, Tp) (l + - g,,{E,) ] F{Ee, Z ^ 1) E, dE, 



(Be 



UEe,Tp) ( 1 + - gr,{E,) + - UEe) ] F{Ee, Z=l)E, dE, 



(1-18) 



For the calculation of the proton recoil asymmetry C, first, we have to integrate over the azimuthal angles. After the 
integration we arrive at the expression 



d'^Xp- {E^, kp,dp, 9ep, P) 
dEgdkpd cos 6pd cos 9ep 



, qN2N I , fcg + fcefcpC0s6>ep , 1 1 

^^ + ^^ ^^^y-"^^^^ eIe + I7TT3A^ 



x\Eo-^ + {X^ + 2{k+1)X 



kl + fcgfcpcos 
E 



- A^-2(k+1)A 



-|-Pcos6'p 



2/1^ \ cos dep r)/ COS 6(,p + fcp 1 1 

^^K^e) 0{Ee) - 



E, 



E 



M 1 -f- 3A2 



Jv^kp cos Q cp 



(2k+ 1) A - A 



fcg ~\~ kp cos Qp 
ELE 



.2(« + i)A h^hv^l^(^kp + fee cos O^p) 



EeE 

Having integrated over cos 0ep we obtain the following result 

d^X0-iEe,kp,ep,P) ' 



1 + - gn{E,) F{E,, Z=l) S{mn - Ep - E^ - E) k^E^kl. (1-19) 



'Pa 

dEedkp 



X 2(£;o£;e - mi){Eo - E,)kp + (1 + 3A^)Se(So - K - K) - (1 - y){Ea - E,){E^ ^ mi ~ fc^)fe 



+ (A^ + 2(k + 1)A) E.iiEa - E^V + ^'p - - (A' - 2{k + 1)A) (£^0 - E,){{Eo - E.f - k'^ ~ ki)kp 

1 1 



-PCOS Or, 



A{E,) {Eo - E,){{Eo - E,y - fc^ - k^) - B{E,) E,{{E,, - E,Y + fc^ - k'^) + — - 



3A2 



X - (2k + 1) A 2(^0 - E,)E,e^ - X {{Eo - E,Y - fc^ H- ki)k; + 1)A {{{Eo - E,Y - k^Y - k^ 



A(l + A) (3(^0 - E,Y -k;- ki){E^ - mi - fc^) + A(l - A)2(So - E,)E,{E'^ - mi - fc^) 



{l + ^g,,{ES)F{E,,Z ^l), 



(1-20) 



where we have kept the terms of order 1/M and the radiative corrections of order a/n. The integration over the 
electron energy Ee and the proton momentum kp gives one 



dX^-{ep,P) _ Gl\Vud\^ r.. a 



dcosi 



(1 + 3A^) 



167r3 l^^ + ^^^ + X7TT3A^ 



X3 + (1 + 3A2) {X4 + Yi) - (1 - A2) (X5 + Y2) 



+ (^A^ + 2{k + 1)A ) - ( A^ - 2(k + 1)A ) X 
1 1 



P cos 0n 



- {Ao + Bo) Xs + A0X0 + A0- Xio -Bo- Xn 

TT TT 



M 1 + 3A2 



(aXi2 -{k + 1)AXi3 - {2k + 1) AXi4 - a (1 + A) (X15 ^3) + A (1 - A) {X^o + n))] }, (1-21) 
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where we have denoted 



X. = 



Xa = 



X. 



Xa = 



/■{fep 



r(E. 



2{Eo - Ee)EekpF{Ee, Z = 1) dEedkp = 0.235044 MeV^ 

/■(£e)max 

/ 2{Eo - Ee)Eekpg„{Ee)F{Ee, Z = I) dEedkp = 3.932276 MeV^ 

(Be).„,„ 

2{EoEe - ml){Ea - Ee)kpF{Ee, Z = 1) dEedkp = 0.225101 MeV^ 
Ee{Ei -ml- kl)F{Ee, Z ^ I) dEedkp = 0.308674 MeV^ 

/(Be).„,„ 
/-(£;e)max 

/ [Eq - Ee){El -ml- kl)kpF{Ee, Z ^ I) dEedkp = 0.112937MeV'^ 



HE. 



riE. 



/(Be 



EeiiEa - Eef + kl - kl)kpF{Ee, Z = I) dEedkp = 0.112937MeV^ 



Xv = 



Xa = 



/.{fcp)max f{E,U^^ 

/ / {Eo-Ee){{Eo-Eef -kl-kl)kpF{Ee. Z= I) dEedkp = -Q.ll2lMyicV'^, 

■Jo "'(-Ee)mi„ 

Ee{{Eo - Eef + kl- kl) F{Ee, Z^l) dEedkp = 
{EQ-Ee){{Eo-Eef -kl-kl 



^(fep)max /-(Be),, 

J{E,)^i 

/■(fep)„ax r(E^)„ 



) + Ee{{Eo - Eef + kl - kl 



0.129702 McV^ 

F{Ee,Z ^ 1) dEedkp = 



= 0.001061 MeV^ 

/•(fcp)„ax /-(Be),, 

Xio - / / 

•/O "'(-Ee)mi 

= -2.157114 MeV^ 



\Eo - Ee){{Eo - Eef -kl- kl){gn{Ee) + fn{Ee)) F{Ee, Z ^ I) dEedkp = 



Xt2 — 



/.(Be)„ax 

/ EeiiEo - Eef + kl- kl) gn{Ee) F{Ee, Z ^ 1) dEedkp = 2.244201 MeV^ 

'(Be).„.„ 
/-(Be)„ax 

/ [{EoEe - ml){{Eo - Eef + kl - kl) - Eo{Eo - Ee){{Eo - Eef -kl- kl) 

Jq JiE^)^i„ 

-{{Eo - Eef -kl + kl)kl] F{Ee,Z = 1) dEedkp = 0.306001 MeV^ 

/.(fcp)max /.(B.)„ax 



-''^13 — 



X, 



X, 



r{E. 

HE. 

r(E 



{{{Eo - Eef - klf - k^)F{Ee,Z = I) dEedkp = -0.119183 MeV^ 



/•(fcp)„ax /•(Be)„ax 

/ / 2{Eo - Ee)Eekl F{Ee, Z ^ I) dEedkp = 0.186818MeY\ 

Jo "'(-Ee)mi„ 

/■(fcp).„ax r{E,)^^^ 

/ / {3{Eo-Eef -kl-kl){E^ -ml- kl)F{Ee,Z = 1) dEedkp ^-0m2705McY^ 

Jo "'(-Be)mi„ 

riE 



Xie — 



Xn — 



/■(fcp).„ax /.(Be)„ax 

/ / 2{Ea- Ee)Ee{E^- ml- kl)F{Ee,Z = 1) dEedkp = 0.291295MeY 

Jo "'(-Ee)mi„ 
/-(fcp).„ax /-(Be)„ax 

/ / -^{{Ea- Ee)^ -kl-kl)F{Ee,Z = I) dEedkp ^-0.211978MeY\ 

Jo "'(-Ee)mi„ 
r(fcp)„ax |.(Be),„a.- 



^18 — 



/.(Be),„ax 

/ me ((i;o - £;e)^ - A:e)-P'(^e, 2- = l)(i^;ecifcp = 0.093562 MeV^ (1-22) 

"'O "'(-Be)mi„ 

The factors X17 and Xig are related to contributions of interactions beyond the SM (see section IIX| . The factors Xj 
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for j = 1, . . . , 16 are calculated for the limits of the integration over Ee and kp equal to 

< kp < (A:p)max = y^^^-m^, 
_GE^_M^+m| _{E„ + kp)^ + ml 



2(So - kp) 



2(-Bo + kp) 



(1-23) 



The factors Yj for j = 1, 2, 3, 4 are defined by the contributions of the 1/M corrections to the limits of the integration 
over Ee and kp. The 1/M corrections to the limits of the integration over E^ are given in Ea. p-12[ ). The limits of the 
integration over kp have no 1 /M corrections 



{kp)ma.^ = \ —{{mn - VTLpY - ml) = J E^ 
V 



(1-24) 



The account for the 1/M corrections to the limits of the integration over E^ can be carried out by the formula 



(-Ee),„i„+A_E_(fcp) 



^{Ee, kp)dEedkp 



(Be),„ax(fcp) 



(Be),„.„(fcp) 



$(£'e, kp)dEedkp 



(kp 



$((Se)max,fcp) AS+(fcp) - $((Se)„,i„, fcp) AS_ (/Cp) 



dkri 



(1-25) 



The factors for j = 1, 2, 3, 4 are defined by 



Yi = 



1 



(^Eq — {Ee)ma.xj (^^e)max ^1 — 



^^((Se)max,^= 1) 



Eq — (£'e)min^ (£'e)min(^l 



(£^0 - fcp)^ 



(£^0 + fcp)2, 

F{{Ee)rain, Z = 1)1 (iJ^ - - fc^) fcp dkp ^ 0.077953 McV^ 



1^2 = 



{kp)„ 



Eo — (£^e)max j ^ + (£'e)^ 



m:f 1 - 



(£^0 + fcp)2 



X F((£;e)max, Z = 1) - ((Sq - (Se)mi„) ~ + (Se)^i„ " (l 



(^^0 - fcp)^ 



X F{{Ee)nun, Z =1) {El ~ml~ kl) kp dkp = -0.000773 MeV^ 



Y. = - 



(So - (i?e)max) ( (i?0 - (i?e)max)' " - (i5,)^^^ + m^) (l - 



X i^((i?e)max, Z ^ I) - (e^ - (i?e)mi„) ( (i?0 - (i?e)mi„)' " fc^ ~ (i?e)fni„ + (l " (^^"""^fe^p ) 

X i^((£;e)min,^ = 1)1 (-Bq - - fc^) ^ -0.134228 MeV^ 



mf 1 - 



(£^e)max( ( £^0 - (-Be)max] + kp - [E^) 

max ' 

X F((£;e),„ax, Z=l)- iEe),r.in[[Eo - (£;e)min)' + fcp " (£^e)Ln + "^e) ( 

X F{{Ee)nun, Z = 1)1 {eI - mj - fc^) dkp = -0.143549 MeV^ 



[Eo + kp)^ 



{Eo - kpf 



(1-26) 



The correctness of our calculation of the electron-proton energy distribution a{Ef,,Tp) we may verify by calculating 
the lifetime of the neutron. Having integrated Ea. p-21| ) over cosSp we obtain the rate of the neutron /3~-decay 



(A^-)sM = (l + SA^)^^'^"^'' 



{Xi + ^ X2 + ^ [Xs + (1 + 3A2) {Xi + Fi) - (1 - A^) {X, + Y2) 



87r3 I TT il/ 1 + 3A2 

+ (a^ + 2{k + 1)a) Xq - (a^ - 2{k + 1)a) Xt\ }, 



(1-27) 



related to the lifetime as (A^- )g-^ = (t„)sm- The numerical factors Xj and Ya; for j = 1, . . . , 7 and fc = 1, 2 define the 
lifetime of the neutron with the account for the radiative {X2 7^ 0) and the 1/M {Xj 7^ for j = 3, . . . , 7 and Yfc 7^ 
for k — 1,2) corrections. 
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Let us consider three possibilities: 1) = for j = 2, . . . , 7 and = for fc = 1, 2 (i.e. without radiative and 
1/M corrections), 2) Xk ^ for fc = 1, 2 and = for j = 3, . . . , 7 and Yfe = for fc = 1, 2 (i.e. with radiative 
but without 1/M corrections) and 3) 7^ for fc = 1,...,7 and 7^ for fc = 1, 2 (with radiative and 1/M 
corrections). For these cases we obtain 1) (t„)sm = 915.3 s, 2) (r„)sM = 881.0s and 3) (t„)sm = 879.6s, respectively. 
These results agree well with the lifetimes of the neutron, adduced in Table I (see section IX)) . 

The correlation of the neutron spin and the proton 3-inomentum is described by the angular distribution (7ll | 



47r 



1 + 2PCcos( 



(1-28) 



where P is a neutron polarisation. The correlation coefficient C is defined by 

1 



C ^ --^{Ao 
2X1 ' 



„ , 1 Xq , a 1 X-i n , ct 1 X-i 1 „ 1 
Bo) + -Ao + Aa Bo + 

"^ 2X1 ' TT 2 Xl ' TT 2 Xi " Ml 



f 3A2 
a 1X2 



A (k + 1 A 

2 Xi ^ ' 2 Xi 



1 



1 



AX. 



TT 2X1 M 1 + 3A2 \2Xi 



+ (1 + 3A2) 



IX 



+ (A2 + 2(k+1)A) 



IXe 
2X1 



-(A2-2(k + 1)A) 



2XJ' 



(1-29) 



where the factor xc = Xg,/2Xi = 0.27491 agrees well with the factor xc = 0.27594, calculated by Gliick [tI]- The 
appearance of the term ^0X9/2X1 is related to the deviation of the Fermi function F{Ee, Z = 1) from unity. 

The contributions of interactions beyond the SM (see Appendix G) changes the electron-proton energy distribution 
a{Ee,Tp) as follows 



a{Ee,Tp) 



1 + 

i\CA\ 



be 
2 



1 



1 



2 1 + 3A2 



\Cy\^-l) + i\Cv\'--l)-i\CA\'-X') 



A^) 



1 



En 



1-A2 M 



{1 + -UE,) C2{Ee,Tp) 



(1-30) 



Multiplying a{Ee,Tp) by the lifetime of the neutron r„ (sec Eq. (|G-lll) ) we obtain 



Tna{Ee,Tp) = (r„)sM|(] 
~i\CA\' ~ \') - {\Ca\' - 



bp-p, bp 



r)]-a^[bc + b 



m, 

Ee I SM 

m 



Cl{Ee, Tp) 

n 



Ep / SM 



ao 
1 



1 1 

2 T 
E, 



3A2 



(|Cyp-l) + (|Cvp-l) 



1-A2 M 



j^l + ^f^{Ep))C2{Ee,Tp)}. (1-31) 



This allows to reduce the analysis of contributions of vector and axial-vector interactions beyond the SM to the 
electron-proton energy distribution to the analysis of these contributions to the axial coupling constant A (see sec- 
tion |VI] and section HXj) . Contributions of interactions beyond the SM to the proton recoil asymmetry C are given in 
section IIXI 
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